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• Abstract 

^ ' We consider the Boltzmann equation for a gas in a horizontal slab, subject to 

a gravitational force. The boundary conditions are of diffusive type, specifying the 
wall temperatures, so that the top temperature is lower than the bottom one (Benard 
^ ^ setup). We consider a 2-dimensional convective stationary solution, which is close 

■ for small Knudsen number to the convective stationary solution of the Oberbeck- 

CN . Boussinesq equations, near above the bifurcation point, and prove its stability under 

2-d small perturbations, for Rayleigh number above and close to the bifurcation point 
and for small Knudsen number. 
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1 Introduction 



We study the small Knudsen number solutions to the Boltzmann equation in a slab with 
diffusive boundary conditions in the presence of a gravitational field. 

In a previous paper [AEMN] we have proved the existence and stability globally in 
time for small Knudsen number of a positive one-dimensional stationary solution to the 
Boltzmann equation, which is close to the hydrodynamic laminar solution of the Oberbeck- 
Boussinesq (0-B) equations. At the hydrodynamical level there is a bifurcation phe- 
nomenon: when the vertical temperature gradient exceeds a certain critical value, the 
laminar one- dimensional solution loses stability and various two- or three-dimensional pat- 
tern flows appear. In particular, it has been proved the existence of a two-dimensional roll 
solution of the 0-B equations close to the bifurcation point. Its stability under suitable 
perturbations has also been proved. 
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In this paper we construct, by means of perturbative arguments (expansion method), 
for small Knudsen number, a positive two dimensional solution to the stationary Boltzmann 
equation, which is close to this roll solution. Moreover, we prove its stability for long times 
under a suitable class of two dimensional initial perturbations. These results are true for 
values of the Rayleigh number above and close to the bifurcation value, provided that the 
force is small enough. To state our result, we need to introduce some notation. 

Consider a gas in a 2-dimensional box of height 2TTd and length 2'n'h, under the action of 
a gravitational force g in the direction z. The upper and lower walls are kept at temperature 
T+ and T_ respectively, with T+ < T_, with no-slip conditions, while periodicity is assumed 
in the horizontal direction. At the kinetic level, the behavior of the gas is given by the 
following Boltzmann equation with boundary conditions diffusive in the z direction and 
periodic in the x direction, written in dimensionless form. 



where 



dF 1 dF 1 dF ^dF 1 

at e ox e oz ov^ e 

F{0, X, z,v) = Fo(x,z,v), {x,z) E {— fin , iJ.n) x {—n,n), v E M^, (1.1) 

F(t, x, =F7r, f ) = Mq;(t>) / \wz\F{t, x,^n,w)dw, t > 0, ^ 0, x G [— /xvr, /itt], 

Fo>0, M_ = ^e-^, M,(.) = ^-^^-L_^e-^, 

Qif,9){z,v,t) = ^ I dv^ I dujB{u,\v-v^\){flg' + f'g',,~ f^g-gj}. 
Jr^ Js2 

Here h',h'^,h,h^: stand for h{x, z,v' ,t), h{x, z,vl,t), h{x, z,v,t), h{x, z,v^:,t) respectively, 
5*2 = {a; e |ci;^ = 1}, B is the differential cross section 2B{lj, V) = \V ■ u!\ correspond- 
ing to hard spheres, and v, and v' ,v'^ are precollisional and postcoUisional velocities or 
conversely. The boundary conditions are such that the condition of impermeability of the 
walls, 

dvFvJv = 0, (1.2) 



is satisfied. The solution depends on the parameter e = 2\/6Kn/\/5TT, where Kn is the 
Knudsen number given in terms of io, the mean free path of the gas in equilibrium at 
temperature T_ and density p. We have also put Ma = €\/6/ -\/5, where Ma is the Mach 

number. With this choice, the Rayleigh number ([Soj) is Ra = ^ ^ , independent 

IT 

of e. We will fix through the paper the parameter G such that G < Gq with Gq suitably 
217 d 

small. Fix h = where ac is the critical wave number for the first bifurcation. The 

linear analysis of the 0-B equations with rigid-rigid boundary conditions, (13.11) below. 
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describing the behavior of the fluid at hydrodynamic level, gives a critical value ac and 
a corresponding critical Rayleigh number |DRj . The parameter A will be chosen 
in an interval [Ac, (1 + 5) Ac], for 6 small and with Ac determined by the condition that 

Rttc = -^^1 — is the critical value. With these choices, at the hydrodynamical level, 
vr 

two roll solutions will appear at the bifurcation point, consisting of one roll, rotating 
clockwise and anticlockwise respectively. These solutions are constructed perturbatively 
for 6 small in a rigorous way [lu] and their local non-linear stability has been proved for 
small initial perturbations with the same period of the roll solution pu] . 
The clockwise solution hg is then of the form 

K = he + 5 Kon + 0{S^), (1.3) 

where hi is the laminar solution and hcon is the eigenfunction corresponding to the least 

eigenvalue of the linearized Boussinesq problem around the laminar solution (see Section 

4 for the precise definition). 

In this paper, we construct a stationary solution Fg of the Boltzmann equation, which is 

close for e small to the hydrodynamical solution (say, the clockwise one) in the sense that it 

1 _^ 

can be written as a truncated expansion in e, Fg = M + efs + 0(e'^) with M = - — 2 
and 

/ If P — " 
fs = Myps + Us ■ V + Ts — 

where p^, Ug-, Ts are expressed in terms of hs. Moreover, we prove the kinetic non linear 
stability of Fg under suitable initial perturbations. 

We study the Boltzmann equation for the perturbation $ = M^^{F — Fs) with the 
initial datum 

5 

^o{x, z,v) = J2 ^""^^"^ (0' a;, z, v) + e% (1.4) 

n=l 

where ^^^\0,x, z,v) is the n-th term of the expansion introduced in the next paragraph, 
computed at time t = 0, and the e-dependent is arbitrary but for having total mass 
J dvdxdzM{v)p5{x, z,v) = and satisfying (3.8). 

We write also the time dependent solution in terms of a truncated expansion in e 

5 

^{t,x, z,v) ='^£"'^^''\t,x, z,v) + eR{t,x, z,v), {x,z)eil^, (1.5) 

n=l 

where = [— /ivr, yuvr] x [—it, tt]. The first term of the expansion in e is 



where the fields p^(t, x, z),u^(t, x, z), 9^(t, x, z) are solutions of the hydrodynamic equations 
for the perturbation, with initial datum (-uj, ^q). The initial data are chosen as follows: let 
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{uq,6q) be an initial perturbations of the convective solution {us,Os) sufficiently small to 
ensure that the solution {u(t,x, z),9(t,x, z)) = {us{x,z) + u^{t,x, z),9s{x, z) + 9^{t,x, z)) 
of the initial boundary value problem for the 0-B (13. ip equations exists globally in time 
and converges to {us, 9s) as t — > +00. 

The next orders involve kinetic corrections in the bulk as well as boundary layer cor- 
rections. In fact, at next orders, the standard Hilbert expansion bulk terms do not satisfy 
the diffusive boundary conditions and boundary layer correction terms are to be included 
to restore the boundary conditions. They are computed solving suitable Milne problems 
in the presence of a force F and a source term S. Under suitable assumptions this prob- 
lem has ^-smooth solutions when F is a potential force decaying fast enough at infinity 
|CME] . We give in Section 3 a procedure to compute the $„'s in the time dependent case 
and show that they have good enough properties as consequence of the smoothness of the 
hydrodynamic solution. In particular, they inherit the smallness and decay properties of 
the hydrodynamic solution, such as the exponential decay in time. The main difficulty is 
then the control of the remainder R asymptotically in time. The equation for the remain- 
der i? is a weakly non-linear Boltzmann equation with a source B, generated by the terms 
of the expansion. Since it is weakly non-linear, it is enough to get good estimates for the 
associated linear problem, which is of the form 



where L is the linearized Boltzmann operator, defined in Section 2, and J{^h,R) is a 
linear operator depending on the perturbed stationary solution and on the first terms of 
the expansion. The operator L is non positive on L2(M^, Mdv), but has a non trivial null 
space Kern(L). The linear operator J{^h,R) is at least of order e but it is the main 
contribution for R G Kern(L). The control of the component of R in the orthogonal to 
Kern(L) is given by the well known spectral inequality for L (see e.g. |Maj ) . 



where (/, Lf) is the scalar product in L2(M'^, Mdv), P is the projector on Kern(L), z/(f ) = 



Using this (and ignoring the boundary contributions) one gets the differential inequality 



with II ■ ||2,2 the norm in L2(fi^ x M.^,Mdxdzdv). This produces bounds growing expo- 
nentially in time. To avoid this we use a spectral inequality for the operator 
Lj{f) = Lf + eJ{^Hi Pf)- The inclusion of the second term in the new operator simplifies 



di 



-1 O O r\ -, 

R + -K^^ + ^^Q-^R) - M-'G—{MR) = - 



LR+-J{<^H,R) + B, 



{f,Lf)>C{{l-P)f,u{l-P)f) 



(1.6) 
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all the arguments in the control of the remainder. The price to be paid is that Lj is not 
self-adjoint and its null space is more complicated. This inequality is an important new 
ingredient for proving stability results in the Boltzmann equation framework. In [ANlj 
there was no need to use this inequality to get the global stability result because in that 
paper it was possible to take advantage from the fact that the constant C in the previous 
inequality can be taken small, while the parameter controlling the bifurcation can grow 
beyond the bifurcation value. In the Rayleigh-Benard setting, there is not such freedom. 

To study the hydrodynamic part of R, PR, we use a duality argument which was 
introduced by N. Maslova [Ma] . To illustrate how it works, let us consider the typical 
equation one has to study in the stationary case 

vV.R- eG^^iMR) = -LjR + q 

with prescribed incoming data ai z = ivr, periodic in x. One also considers the dual 
equation 

z; ■ - £G^^(M0) = h}^ + h, 

with vanishing incoming data at z = ±7r, periodic in x. By taking the inner product of 
the first by (p and the second by R and summing, one gets an estimate of \\h\\ in terms of q 
and with suitably small coefficients. For that, we need G small, but we can anyway have 
large Rayleigh number by increasing A. Then one chooses h = PR and thus shows that 
1 1 Pi? 1 1 can be bounded in terms of ||0|| and The equation for is studied by means of 
Fourier analysis, which provides a good control on but for the 0-moment. The estimate 
of the 0-moment is based on a direct approach, by means of rather explicit calculations of 
the first few moments of the solution using o.d.e. analysis, for the one-dimensional case. 
This allows a reasonably simple analysis when the incoming data are prescribed. Dealing 
with diffuse reflection requires more technical steps. This method has been exploited in 
[ANlj and has been extended in |AEMNj to the one-dimensional Benard problem, which is 
more difficult to deal with because of the presence of the force and the diffusive boundary 
conditions. In this paper, we extend the one-dimensional analysis to the two-dimensional 
case by using that, both in the conductive and convective case, for 7?. < 7lc{l + S) and 
S small, the contributions due to the inhomogeneity in the x variable are small and can 
be included perturbatively. A by-product of this analysis is the extension of the result in 
[AEMN] to two-dimensional initial perturbations. 

Finally, the control of the nonlinearity requires L°°-estimates (in space), which are 
more intricate by the presence of the force. We use techniques based on the study of the 
characteristics, which are no more straight lines because of the presence of the force. This 
is controlled by looking at the characteristics in the phase space like in |EML] and [AEMN] . 

The main result of this paper is summarized in the following theorem. 

Theorem 1.1. Assume that the gravitational force is such that G < Gq, with Go small 
enough, and X G [Ac, (1 + 5)Ac]. Then, there are 6o and Eq small enough such that for 
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6 < 6o, there exists a positive steady solution Fs to the Boltzmann equation, such that for 

II M-'[Fs~{M + 6fs)] h2<ce^ 

Assume also that the initial perturbation matches the expansion up to orders^ as detailed in 
Section 3 below, and is small as detailed in Section 4. For such perturbations the stationary 
solution is stable uniformly in e for e < Sq. 

Here, stable means that the perturbation vanishes asymptotically in time. This is a 
consequence of the inequality 

/ dt dxdz / dv\^{t )| M{v) < oo , (1.7) 

Jo Jn^ Jr^ 

which is proved in Section 4, and of the regularity of the solution which follows from our 
construction. 

The positivity of the stationary solution is obtained by using the methods in |AN4j . We 
remark that the method presented here for proving stability strongly relies on the fact that 
the problem we are dealing with has suitable stability properties at the fluid dynamic level, 
which we show to be preserved in the kinetic setup by means of a perturbative analysis 
starting from an Hilbert-type asymptotic expansion plus boundary layer corrections. The 
preservation of the fluid dynamic stability at kinetic level also occurs in the Taylor-Couette 
case discussed in |AN1] . where the bifurcation phenomenon also arises. 

The paper is organized as follows: in Section 2 we construct the stationary solution 
as a Hilbert asymptotic series in e. For sake of shortness, we choose not to give here the 
construction of the terms of the expansion, and refer for that to Section 3, where we show 
explicitly the analogous construction in the time-dependent case. We do complete the proof 
of the existence of the solution in the stationary case in Section 2, by proving the main 
theorem on the remainder. In Section 4 we deal with the remainder in the time-dependent 
case and prove the stability result. 



2 Stationary case 

In this section the stationary case is treated. Here, (and also in Section 4), for sake of 
simplicity, we consider a square box [— vr, vr]^ instead of a rectangular box [—fin, fin] x [—n, n] 
and hence the x-derivative in the equation will have a factor ft in front. 

We write Fg = M(l + in terms of a truncated expansion in the Knudsen number e 
plus a rest term: 

5 

^l{x,z,v) = ^e^^i^\x,z,v) +eRs,eix,z,v). 
1 

The expansion will not be given explicitly here since all the ideas and details are in the 
previous paper |AEMNj . We only remark that the construction of the 's relies on the 
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solution to a Milne problem with external force given in [CMEj . Moreover, we will give 
in the next sections more details on the analogous expansion in the time dependent case. 
In this section we study the equation for the remainder, beginning with some results on 
linear existence together with corresponding a priori estimates. The section ends with an 
existence proof for the nonlinear stationary rest term. 

We denote hy H = L\^{^^) the Hilbert space of the measurable functions on with 
inner product (■,■) = (■, M- )2, where ( ■ , ■ )2 is the standard inner product and || ■ II2 
the standard L^-norm, while || ■ || is the norm corresponding to ( ■ , ■ ). 

The linearized Boltzmann operator is defined, for any / in a dense subset of H as: 

Lf = 2M-^Q{M,Mf). (2.1) 

It is well know that it can be decomposed as L = — z/J + K, where / is the identity, K 
a compact operator and defined in the Introduction, is a smooth function satisfying 
the estimates z/o(l + \v\) < < z/i(l + |f |) for some positive z/q and ui. The operator L 
is a non positive self-adjoint operator with domain = {f & H \ ||z^2/|| < +cx)}. 

The functions i^o = 1, V^i = = v^, i)2 = i^y = Vy, i^s = i^z = v^, i)^ = - 3) fo™ 

an orthonormal basis for the kernel of L in if, Kern(L). For any function in H introduce 
the orthogonal splitting f = f\\ + f± Pf + {I — P)f, where /y is called fluid dynamic 
part and is given by 

4 

fll{x,z,v) = ^fjipj, fj ■= {f,iljj), j=,0, ...,4, 

j=0 

while the non hydrodynamic part f± satisfies ijjj) = 0, j =, 0, . . . , 4. P denotes the 
projector from H on Kern(L). Note that the range of L is (/ — P)H. We will use the 
same symbol for the projections also when dealing with functions depending on x, z, t. We 
remind that the operator L satisfies the spectral inequality (II. 6p . 

For 1 < g < +00, let L'^ be the Banach space of the measurable functions from [— tt, tt]^ 
in H, identified with the space of measurable functions from with norm 

2X 1 



that is 



/||,,2= I / dvM{v) { / dxdz\f{x,z,v)\' 

'[-7r,7r]2 



=■■= {/ : [-TT, tt]^ X ^ M I || / ||,,2< +00} 



Moreover ( ■ , ■ )2,2 is the corresponding inner product for q = 2. 

We also need a function space for the boundary functions. We denote by the sets 
V = {vx,Vy,Vz) such that Vz ^ 0. We consider the functions on [— vr, vr] x {— vr} x M.^ U 
[— vr, tt] X {tt} X M.^ and define the norm 




/L,2,~=sup / dv\vz\M{v)[ / dx\f{x,TTr,v)\'i 
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The Banach space L"*" is the set of such functions with finite || • ||2,2,~ norm. The ingoing 
and outgoing trace operators 7^ are defined by 



7^/ 



f\z=n, if W G 



d3 



The function space where the stationary solution will be constructed is the space 

W^'- := {/ : [-71, 7i]^ xR^ ^R\u^2f e ly-^Df G 7+/ G L+}. 

Note that the norm || ■ ||2,2,~ is defined only for incoming velocities. In the sequel, with 
an abuse of notation we will denote by || 7"/ ||2,2,~ the || ■ ||2,2,~-norm of 5*7"/, where S 
is the reflection of the z component of the velocity. 

We do not explain here how to construct the terms of the expansion We simply 
state a theorem about their properties. We assume that the Rayleigh number Ra is in 
{Rac, Rttc + S) with 6 > and sufficiently small and will consider, for sake of definiteness, 
the clockwise convective solution corresponding to it. 

(n) 

Theorem 2.1. The functions $s , n = 1, . . . , 5 and ipn,e can be determined so as to satisfy 
the boundary conditions 

+tpn,e{^^ =F7r, W), t>0, V^^O, 

and the normalization condition J^^^^_^ ^^^2 dvdxdz^^"'^ = 0, so that the asymptotic expan- 
sion in e for the stationary problem U.l\) . truncated to the order 5 is given by 

F^'''P\x,z,v) = M{v) i^ + ^t'^^'^\x,z,v)^ . 

The functions ^'•"^ 's satisfy the conditions 

II ||2,2< 00, II ||oo,2< 00 , n=l,...,5. 

Moreover the ^^"'^ 's differ from those of the laminar solution by 0{6). 

The functions ilJn,e o,i"c such that ||'?/'n,e||(j,2,~j q = 2,oo are exponentially small as e —>■ and 

/jg3 dvVzM{v)ipn,£ = 0. Finally, there exists a stationary solution to U.l\) in the form 

F^ = F('-p)+eRs,e. 
The remainder Rs,s, simply denoted by R solves the boundary value problem 

Ixv^-^R + Vz-^R - sGM-^ ^^^^^ = -LR + Y^'~'Ji^^'\R) + JiR,R) + A (2-2) 
ox oz ovz e ^-^ 

1=1 
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R{x, =F7r, v) = ^j/^^^ [ \w,\M{w)(R{x,T'^,w) + -iJe{x,TT^,w))dw--iJe{x,T'^,v), 
M{v) J^^^o V e J e 

for f 2 ^ and x G [— vr, it] (2.3) 
2 

where ^t/'^ = — X]n=i ^"^V'n.s? fi') = J^^^^^^ ^d) ^'^^ A is a smooth function hounded 
in II ■ ||g^2; = 2, oo, sfic/i that I dvM{v)A = 0. 

We do not give here the exphcit expression of A for which we refer to |AEMNj . It has 
to be considered as a known term in the rest of this section as well as the ^'^■'^'s. The 
terms of the expansions contribute to A together with their space and velocity derivatives. 
It suffices to know that it has finite || ■ ||g,2-iiorm. Note that M and M_ differ just by 
the normalization and that the ipn,e^ are known functions exponentially small in e due 
to the boundary layer corrections in the expansion. With these boundary conditions the 
remainder satisfies the impermeability condition 



dvv.MR = (2.4) 



at the walls. 



Now we start the analysis of the equation for the remainder R. The first important 
result is a spectral property for the operator Lj below: fix x, z and define, for each f & H 

Ljf = Lf + eNPf , (2.5) 

with 

5 

iV = J(g, ■), q = J2e--'^^:\ 

n=\ 

Note that in the time-dependent case the form of q will be slightly different but without 
affecting the argument below. 

The operator Lj is not symmetric in H. We denote by L} its adjoint. To characterize 
the kernel of Lj, Kern(Lj), note that the functions 

^\)j = ^\)^-eL-^N^\)^, (2.6) 

are in Kern(Lj), where, for G (/ — P^H , L^^f denotes the unique solution of Lg = f 
orthogonal to Kern(L). In fact, Ljipj = Lipj - eNPipj + eNP^jj - e^N P[L^'^N Pi)j] = 0, 
because ipj are in Kern(L) and the last term is zero since P kills the terms in the range of 
L-\ 

It is easy to check, by using (11.61) . that, at least for e sufficiently small, they actually span 
Kern(Lj). Indeed, suppose that there is g E Kern(Lj) with ||5f_L|| = 1 such that {g, ipj) = 0. 
This implies = eY.^j^oi9±, L~^N^|Jj)^pj. Moreover, = ig,Ljg) = {g±, Lg±)+e{g±, Ngu). 

Therefore by (USD, C < -{g±,Lg^) = e^Y.U^9±.Nij^){g^,L-^Ni;,) < ae\ for some 
positive a independent of e. This is in contradiction with C > for e sufficiently small. 
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Let Pj denote the orthogonal projection onto Kern(Lj). Since the range of Lj is 
(J - P)H, it follows that Kern(L}) = PH = Kern(L). 

Consider the space J\f generated by {ipo, . . . , ip4, L'^Nil^o, . . . , L'^Nip^}. We decompose 
this space into Kern(Lj) and its orthogonal complement in Af, Li. 

Moreover, (/ — P)H can be decomposed in the span of {Nijji, . . . , Nil){\ and its orthog- 
onal complement, Lq. It is easy to check that 

H = Kern(Lj) © Li © L^. 

Indeed it is enough to show that Nipi^C. = 0,...,4 is a combination of L^^Nipj^s and 

4 

Ni)(, = ^at,jL''^Ni)j+u, {NiJk,u) = 0, e,k = 0, . . . ,4. 

j=0 

We take the inner product with Nipj. 

4 
j=0 

The matrix with elements {Nipj, L~^Nipk) is non singular, hence the aej are uniquely 
determined, and -u as a consequence. 

Proposition 2.2 (Spectral gap property of Lj). There is Eq > such that, for < e < sq, 
there is c independent of e, {x,z) and 6, for which the following inequalities hold: 

-{Lj^, ^) > c{u{I - Pj)v9, (/ - PM, (2.7) 
-(L>, ^) > c(z/(/ - P)v, (/ - P)v) (2.8) 

Proof . First take (x, z) and the Rayleigh number fixed. It is then enough to consider the 
set H of all functions ip = {I — Pj)ip = aipo + b^i, where ipi e Li and = (J — P)v?o 
G Zq, and with v^ipo and of norm one and + 6^ = 1. First notice that {v(p,Lp) is 
uniformly bounded in Ti. It remains to show that the left-hand side has a positive bound 
from below. But can be decomposed as a sum of an element ipn in Kern(L) with norm 
of order e, and an element L~^Nipi2 in the span of {L'^NipQ, L'^Nip^}. Then 

-{Ljip, V?) = - (aLipo + bN(pi2 + beNipn, aipo + bipu + bL'^N(pi2^ 

The first equality follows from the fact that ipu is in Kern(L) and is orthogonal to the 
range of L; the second equality is due to the fact that {fo,Nipi2) = = {ipQ,Nipii) by 
the definition of Lq; the bound (11. 6p has been used to obtain the inequahty. Since the last 
term is of order e^, it follows that for e > and small enough 

C 

-(Lj(^,^) > -{a\u^±,^^) +b\L-^Ny^i2,L-^N^i2)) > c 
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for some c > 0. The first inequality in the proposition then follows, since the constant 
depends continuously on (x, z) and the Rayleigh number is in a compact set. The second 
inequality is obtained by similar arguments. □ 

We only consider Rayleigh numbers in a suitably small neighbourhood of the first 
bifurcation point, and take G sufficiently small as specified in the sequel. Constants which, 
independently of the parameter e, can be made sufficiently small for the purposes of the 
proofs, will generically be denoted rj. We will first give estimates in the linear case. The 
argument is inspired by the approach in [Ma] and heavily relies on the study of the space 
Fourier transform of R. We use the following definition of Fourier transform: for any ^ E 1^ 



When we want to specify that the Fourier transform is taken with respect to the vari- 
able X we write JF^. In the sequel, if / is a function of {x,z) G [— vr, tt]^, f{^x,Cz) = 
{J-'x^zf){^x, ^z)- Finally, if / is a function of x, z and v, we define < / > as the zero order 
Fourier coefficient of /, that is 



</>:=- — — / f{x,z,v)dxdz, a.a.veM. (2.9) 

(27r)^ y[_,r,7r]2 

and /:=/-</ >. 

An important tool in the analysis is the Green inequality, which we will use extensively 
in the rest of the paper in various different situations. 
Consider the linear boundary value problem 

A..§^+".§^-.GM-'^ = iw + ,. (2,10) 

OX oz ovz e 

with j Mgdv = and prescribed incoming data 

f{x,±7r,v)=p{x,±7i,v) v,^0 (2.11) 

Due to the presence of the force G, the argument requires some care. We introduce 
k{z) =e^^(^+''\ Then we muhiply (2.10) by 2/Mk, integrate over [-7r,7r]2 x and 
integrate by parts to get the Green identity 

II l^h'f Il2,2,~ -^('^/' Ljfh,2 = {fig, f)2,2 + ||«:^P||2,2,~- (2-12) 

Apply the spectral inequality (2.7) to obtain, for small 1] 

II ^h-f\\h,^+^^ II K'^iyHl-Pj)f\\l, 

< c{e II k'^u-^I - Pj)g +77 || n^Pjf +^ || k^Pjq +\\k^p\\1^^^{2.1?,) 
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Finally, since 1 < k(z) < e^^'^*^, we conclude that, for some constant C it results: 



7-/llL,~+7^ II i^Hi-Pj)f\\l 



2e 

< C(^e II u-^I-Pj)g +r] \\ Pjf +^ || Pj9 11^2 +IHl2,2,~)- (2-14) 

Inequality fl2.14p and its variations will be referred to in the rest of the paper as Green 
inequality and used extensively, with L} replacing sometimes Lj. In the last case the 
projector in the r.h.s. has to be replaced by PJ = P. 

Lemma 2.3. Let (p{x,z,v) be solution to 

periodic in x of period 2t{, and with zero ingoing boundary values at z = —tc, tt. Then, for 
some small rj 

II z/5(/ - P)ip ||2,2 < C(^e II - P)g ||2,2 + II Pg h,2 +r]e \\< Pip >||2 ),(2.16) 

II ^ l|2,2 < C( II Z/-5(/ _ p)g 112^2 +^ || Pg h,2 +V \\< >||2 ).(2.17) 

Proof of Lemma 12.31 The method from |Ma] (a variant of |Maj Sen 5.3) can be adapted 
to the present setting to obtain the existence of a solution to (12.151) . if one includes the 
above spectral estimate for L}, and the new characteristics curves due to the force term. 
We write f l2.15p in Fourier variables. For ,^ G — {(0, 0)}, 

i.v^ + izv.)^ = ^Lyp + sGM-'^^^^ +g- {-l)^^\v,\r. (2.18) 



Here 

J^x^i^x^T^^v) for Vz > 0, 
, ^x<^{^x, -TT, v) for < 0. 

Introduce 



(2.19) 



Z = e-'Lyf + g- (-1)«^ \v,\r, Z' = e-'Lyf + g, U={ii- v^)-\ 

Let X be the indicatrix function of the set {t> G | | ■ \< a;|,^|}, for some positive a 
to be chosen later. Let C,s{v) = (1+ \v \y . For ^ ^ (0, 0) 

4 

< c II C-sX$ II Yl II ^^^^i II- II ^-'^'^ 

j=0 
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We use this estimate with the following choice of a, a —\\ C-s^ ||~^|| C-sZ' \\. 

We also introduce an indicatrix function Xi with a = \/Si. We fix Si so that c\/Si << 1. 
Then we find from the above estimate that the P-part of the right-hand side, || P(xi$) ||, 
can be absorbed by || P(xi$) || in the left-hand side. The estimates hold in the same way 
when xi is suitably smoothed around For the remaining x'^x5$ = (1 — x)(l — xi)$ 

we shall use that $ = -UZ. Then 



II pixix'^) r < c II Cs+2x'ix'u Hi c-sZ' r ^ 



+ £6-161 



< 



II c-sZ' r ^ 



II 



+ eG\e\, 



with 



j=0 



6> = E / ^^^^^"^^ 



d{M(f) 



dv 



ij,X'ix\^-UZ)Mdv] . 



We replace a by || C-s^ \\ ^\\ C-sZ' \\ in the denominator. That gives 



p$ ||^<c(||C-.$ IIIIC-.^' 11 + 



\Vz\r 



+ <^i II C-s{I-P)^ \?)+eG\Q\. 



Hence, 



P$ 



Vz\r 



< c((|| II + II C-.(/ - P)^ II) II C-sZ' II 

+5i ||C-.(/-P)*ll')+£G|e|. 



< K II II' + J^|2 + II - nil C-^' 



Consequently, 

II P$ 11^ 



We next discuss the term eGjOj. The first integral can be bounded by e times an integral 
of a product of M, 1 + a polynomial in w, | <^ | and tj oi U^. So this factor is bounded 
by ec II $ ||. And so. 



+ 11 C^s{I-P)<^r)+eG\Q\. 



P$ 



^< c( II C-sZ' 



II {i-p)^ f)- 



Therefore for ^ (0,0), 

\P^^{^,v)Mdv 



<c(l iic-.(^)^(e,-) 11^ + 11 (/-me, 

, II 



+ 11 ^-^^(e,-) ir). 



(2.20) 
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We remember that the zero Fourier mode of P(/9 is zero by definition. Hence, taking e small 
enough and summing over all 7^ ^ ^ Z"^-, implies by Parseval and the spectral estimate 



on L}, that 



+ / i^~^g^ix,z,v)Mdvdxdz+ \\ 7"v9 \\l^2,^ (2.21) 



Inequality fl2Aill implies (l216ll . Replacing (12A4D in flOTll gives (12T711 . and this concludes 
the proof of the Lemma 12.31 □ 

Remark. The statement of Lemma 12.31 still holds if we replace the operator L*j with 
the operator Lj and the operator P with Pj, with some minor modification. The main 
change is due to the fact that one has to use the basis functions in Kern(Lj) namely the 
'i/'/s, instead of the ip/s. They depend on {x,z), therefore we fix a point {xq, zq) and use 
the V^j's computed at this point and the corresponding projector Pjq. Then the argument 
of the proof can be repeated word by word. At the end, we replace Pjo with Pj since 
Pj-Pjo = Oie). 

Put H{R) = Xln=i ^"^^'^('^s"^ -R) decompose H in accordance with the operator Lj. 
Set Hi{-) = H{ ■ ) — J(g, P ■). We notice that Hi{-) is of order zero in e and only depends 
on the nonhydrodynamic projection {I — P). 

To pass from the linear results to the non linear case we will use an iteration procedure 
that will lead to Theorem 12. Ql below. To separate the difficulties coming from the non linear 
term and from the boundary conditions, we split the remainder R in two parts, Ri and R2, 
solutions of two different equations. In the equation for Ri the boundary conditions are of 
given indata type and the nonhydrodynamic known term is included, while in the equation 
for i?2 the boundary conditions are of diffusive type and the known term is absent. The 
equation for R2 will be given later (see eq. (12.241) ). We start with a discussion of the 
equation for Ri, 

/. ...^ + - ^GM-^^-^ = hjR, + H,iR,) + (2.22) 

Ri{x, =F7r, v) = -^^(x, =F7r, v), ^ 0. 

Here Ri is periodic in x of period 27r, and Lj = L( ■ ) + eJ{q, P ■ ) has been introduced 
earlier. An existence proof for this problem can be obtained similarly to that for above. 
The nonhydrodynamic part of Ri is estimated similarly to the corresponding proof for 
Lemma 12. 3^ 

^11 1'Ri Il2,2,~ +^ II i^Hi-Pj)Ri Wh 

< C( II u-^I - Pj)g \\l2 +^ II PjR, \\l2 +^ II PJ9 Wh +^ II ^ Il2,2,~ ), 
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for small rj > 0. Here we have used the fact that 

\{Rl,H,{R^)),^2\<C{\\ u^^{I-Pj)R^ \\ PjR, 

A priori bounds for PjRi will be based on dual techniques involving the problem fl2.15l) . 
Consider first the problem for Ri without the term Hi{Ri). It holds 



Lemma 2.4. Set h := PjRi. Then there is 60 > 0, such that for < 6 < 60, 

1 2 1 

Pj9 2,2 +-i 

Proof of Lemma 12.41 The function Ri is 27r-periodic in x, and here solution to 



h \\l2< C{\\ u-HI - Pj)g +^ II Pjg II ^ III2. 



dRi dRi id(MRi) 1 „ ,^ 

ox OZ OVz s 

=F7r,t;) = -iv^(a;, =F7r,t;), v^'^Q. 
Let V? be a 27r-periodic function in x, solution to 

ox OZ OVz £ 

with zero ingoing boundary values at 2; = — vr, vr. 

We consider the equation for Ri multiplied by Mnip and the one for (f multiplied by 
MkRi and add them. After integrating on [— 7r,7r]^ x and integrating by parts, we 
obtain, 

[ dxdzdvMni^ivzRi^) - eG^^^^^\ = [ dxdzdvMK\-((I - P)<^Lj(I - Pj)Ri) 
J \dz J J le 

+^(/ - Pj)RiLj* (/ - P)if + KQif + nhRi 
Using again the bound 1 < k < e^'^'^^ and the assumption h = PjRi, this gives 
II h \\l,< ^ II j-Ri III2,. +^ II 7> ||2,2,~ +^ II - Pj)Ri Wh 

II '^hi - p)v Wh +^ II '^-hi - Pj)9 \\h +^ II - PW IIL 

+^ II III2 II PJ9 \\h +^ II P^ Wh +^ II \\h, 
for arbitrary positive constants Kj, j = 1, ...,4. It then follows that 

II h \\l,< c[(^ + ^) II ^ |||2,~ H^K, + K, + K,e) \\ z/~l(/ - Pj)g 



, I 1 1 1 ^1 m 7, l|2 ^3 



+ (r/iiri + 7717^3) II h 11^2 +r^(-|- + -i- + -i- + ||<Pv^>||^ 
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For an estimate of the final < Ptp >-term when G is small and the Rayleigh number for 
the rolls lies in a neighbourhood of the bifurcation point, we may apply an exact, direct 
approach based on ordinary differential equations. Namely, < ip >x (■) := / ip{x,-)dx 
satisfies a 1-d stationary laminar given indata problem, similar to eq. (3.5) in |AEMN] . The 
present case is different from the one in |AEMNj because there will be new contributions 
coming from the terms in the expansion depending on x, for example the term q in the 
J-term. Since the first order of the e expansion is of order 5, the same is true for the higher 
order (in e) terms. The new contributions can all be considered as deviations from their 
x-independent values at the bifurcation point, hence are of order 5. We include them in 
the estimate as || Lp \\2,2i V small, and obtain 

||< >h < c ||< Pv >xh,2< c(||< Pip >xh,2 +e II if 112,2) 
c c 

< -\\<h >xh,2 +V II (f ||2,2< - II h ||2,2 +V II (p ||2,2 • 

For details cf Lemma 3.5 in |AEMNj . 

Choosing e « 1, then Ki and K3 {Tesp.K2) of order e^^ (resp. s^"^), rji of order e, and 
using Lemma 12. 3^ the inequality of Lemma 12.41 follows. □ 

Remark. From here on, small factors rj in the estimates will depend also on the small 6o- 

In the following lemma we get the final estimates for Ri. 
Lemma 2.5. If Ri is a solution to the system ( \2.22i\ . then 

II I^^Rl \\2,2 < II Z^~^(/ - Pj)g \\2,2 +^ II Pj9 \\2,2 \\ ^ ||2,2,~ ), 

II l^^Rl ||oo,2 < c(^^ II Z/"^(/ - Pj)g ||2,2 +^ II Pjg \\2,2 +S \\ l^'^g ||oo,2 \\ ^ 

Proof of Lemma 12.51 Consider first the solution to (12.221) with Hi = 0. It satisfies 

II l~Rl ||2,2,~ +£~^ II - Pj)Rl ||2,2< 

c(- II ||2,2,~ II - Pj)g ||2,2 +?7l II PjRl ||2,2 +— II Pjg ||2,2 ), 

for any r/i > 0. By Lemma 2.4 and some additional computations using the solution 
formula, 

1 /I 1 _i , 

II V^Rl ||oo,2< C(^- II V'^Ri ||2,2 +e II V ^g ||oo,2 + II 7^1 l|2,2,- 

It is easy to see that adding the term ifi(i?i) does not change the above results. That 
proves the lemma. □ 
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Now we study R2 , the other part of the remainder. Denote by 
f-{x,T'^,v) = ^j/^.^ / (Ri{x,Ttt,w) + -ilj{x,TTr,w))\w^\Mdw, v^^O, 
the incoming data for R2 which is solution to 

/...^ + - eGM-'^^ = hjR2 + H,m, (2.24) 

ox oz ovz e 

R2{x,TTT,v) = f~{x,T'^,v) + ^j/^? / R2{x,T'^,w)\wz\Mdw, Vz^O. 



M{v) 

Existence and uniqueness for (12.241) follow as in the laminar case, cf |AEMN] . The following 
a priori estimates hold for R2. By Green's formula and the definition of Hi[R2) = Hi{{I — 



II l'R2 |||2,~ +^ II l^Hl - Pj)R2 ||2,2<II 1^R2 Il2,2,~ • 

Here we face a problem of diffusive boundary conditions. The ingoing flow is given for cp 
and for Ri but not for R2. The following bound is proved as in equation (4.23) in pML] , 

=F / VzR^{x,±.n ,v)dv < cer] / \vz\MRl{x,±TT,v)dv 

+— / \v,\M{f-{x,±n,v)fdv. 

The computation in |EML] has to be adapted to the fact that the boundary conditions for 
R2 are not purely diffusive but contain the given data /~. From Green's formula we have 
also 

Cl 



e " 



'5(J-Pj)i?2 ||2,2< j {iv„Rlix,-n,v)) - iv„Rl{x,7r,v)))dx. 
Then, using the previous bound we get for any rj > 0, 

— - Pj)R2\\ln<ceriy" [ \vJMRl(x,±7i,v)dxdv + — \\ f- Wj,^ ■ (2.25) 



We need to estimate the terms involving the outgoing parts of R2 in the r.h.s. in terms 
of ||-R2||2,2 and ||/~||2,2,~- This is done separately at ivr. We start with tt. Consider the 
equation for R2, multiply by KMR2 and integrate in velocity over the region Vz > q- Then 
integrate over space, using a smooth cut-off function x{^) which is in a small interval 
close to —IT and 1 close to vr. Finally, integrate over q, ioi qo < q < and go small enough. 
We get 



f dq f dv f dxVzK{-n)MRl{x,'n,v) <'^\\u^{I - Pj)R2 

J qo Jvz>q J — n ^ 

— dq dv dxdzK(z)x' {z)vzMRI + ceG \\ R 



12,2 



■2 II 2,2 
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The term on the Lh.s. equals 

pO p pit 

C|go| II 7 -R2 ||2,2~ + dq dv dxVzMR\{x,T:,v). 

Here in the last term we can replace R2 by the ingoing boundary data so to estimate it as 
dq dv dx[M^ / dw dxWzMR2{x,TT,w) + Mf {n 

QO J q<Vz<0 J —IT Jwz>0 J —IT 

<c{qo)[\\ I-R2 ||2,2~ + || r Il2,2,~], 

with c(go) = o(|go|)- Adding a similar estimate at — vr gives 

II ||2,2~< ^ II uHl-Pj)R2 \\l2 +C II /- 11^2,- +C II PjR2 11^2 • (2-26) 

Replacing in (12.251) we get 

- II u^I-Pj)R2 \\l2<evC II PjR2 IIL +- II /" Il2,2,~ • (2.27) 
e ' ' erj ' ' 

We shall next prove an a priori estimate for the hydrodynamic part of i?2- First we 
consider the 1-d (x-independent) case and then include in the argument the missing terms 
which, as explained before, will be of order 6. The extension of the 1-d results to the 2-d 
case will be based on perturbative arguments in 6. To take into account these terms of 
order 6 we add to the right-hand side an inhomogeneous term gi with J Mgidv = 0, which 
will be of use later on in the proof of Lemma 12. 7[ 

Lemma 2.6. Let R2{z,v) be solution of the 1-d problem and f~ be defined as before, 

v.^ - sGM-^^^^ = hjR2 + + 9u (2.28) 

oz ovz e 

^2(T7r,t;) = / R2{T7r,w)\wz\Mdw + f-{Trr,v), v.^O. 

U 9w = /k3 dvMgi = 0, then it holds that 

II PjR2 \\l2< 4 II f Wl- + II '''^91 \\l ■ 



Proof of Lemma 12.61 Here, Lj is generated by a function q, independent of x, defined by 
q = q + 0{6). In the same way, in the term Hi we retain only terms of order zero in S, 
which are independent of x, and include the remaining terms in gi. In this way we reduce 
the equation to a 1-d equation with a given term. 

Consider equation (I2.28p for i?2 = ^zR2, the Fourier-transform in ^ of i?2- It satisfies 
the equation 

n;,e.i?2 - eGM-'-^{MR2) = e-^L^2 + iWh) - v^ri-l^ + 9i (2.29) 
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with r(f ) now denoting the difference between ingoing and outgoing boundary values, 

r{v) = R2{7c,v)-R2{-7r,v). (2.30) 

For ^2 7^ we use the method of Lemma [2.31 with Lj instead of its adjoint. Consider the 
ingoing boundary values as known, and follow step by step the proof of Lemma 12.31 with 
obvious changes. Let Ljq be the operator defined in (2.5) with q taken at a fixed {xo,zq). 
Define Pjo as the orthogonal projection on Kern(Ljo). We reach the analogue of (]2.20p . 

j \PjoR2{^,,v)\'Mdv < II C.s{v)L^2{^.,-) ir + II iI-Pjo)R2{^z,-) ir 

^llPjoi^^^^ II f + II u-'^im.m,-) f ). (2.31) 

For an estimate of the r-term, we express it with the = term in the Fourier series for 
(1229]), 

v,r{v) = -L^2i0,v) + eGM~'^{MR2iO,v)) + g\iO,v) + lUR2)iO,v). (2.32) 
Inserting this into fl2.3ip . and summing over ^ 0, results in 

j{P^2)\z,v)Mdvdz <c(^^j - Pj)R2)\z,v)Mdvdz 

+ / u^^gl{x, z,v)Mdvdz + e"^ \\ R2 II22 



We are left with the Fourier component PjoR2{^z) for C,z = 0. Estimate separately the 
(/ — P)-component and the P-component of Pjo-R2(0, ■). For (/ — P)PjqR2{0, ■) we obtain 

II (/ - P)Pjoi?2(0, •) ||< C^dl (/ - Pj)P2 ||2,2 + II PjR2 112,2). 

The P-moments are the more involved and will be discussed each separately. We start 
from the f^-moment of R2{0,v). Multiply fl2.28p by M and integrate over z G [— vr, and 
V. Since gio = 0, we have 

v,MR2{z,v)dv = / f'{-7i,v)v,Mdv . 

Jv^>0 

Given two functions h{v) and f{{-),v) we use the notation fh{-) '■= J dvh{v)f{-,v). 
In particular, for h = ipj, j = 0, . . . , 4, we also use the notation fj. We have 



|4..(o)| 



VzMR2{z, v)dvdz 



<c\\r 



To estimate the moments of -R2(0, v) we use the identity 



R2i0,v) = A - 5^ P2(e.,^)(-1)^% (2.33) 
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where A(f) = tt{R2{tt,v) + R2{—tt,v)), which follows from 

R2{n,v) + R2{-7i,v) = ;i V4fe,^^)(e'"^^ +e-*"«0, 

zvr ^ — ' 

by solving for the ^2 = coefficient. 

We write A = 2ttR2{—tt, v) + TTr{v) for f ^ > and A = 27cR2{tt, v) — 7ir{v) for < 0. 

We consider first the ?/'4-moment of i?2(0,f), denoted by R24. We notice that 

4.2^(0) = ^^24(0) j vyAMdv + Ri^Ai^), (2.34) 

where i?^ = (1 — P)R2 and A and -B are nonhydrodynamic solutions to 

L(f^A) = - 5T), L{v^v^B) = v^v^. (2.35) 

Hence, we are left with the control of -R2t)|yi(0)- Fo^' that, we use (12.331) . A multiplication 
of (12.331) with Mv^A followed by a ^-integration, gives 

-R2i;2a(0) = A„2A - y^-R2^;2A(^z)(~l)^''- 

Let us first consider the contribution A„2^. Using the relation A = 27ri?2(±vr, v)^TTr{v) 
for fx ^0, we notice that the first part it is computed in terms of the outgoing fiow and 
then can be bounded in terms of J dvvlAM f~ and ^qv1AR2{-^'k ■,v)Mdv. Then, we 
multiply (12.281) by MxVzA and integrate over f ^ > (similarly at — vr) , to get the bound 

I / vlAR2{±7i,v)Mdv |2< r/ II PjR2 III2 +c(^ || {I - Pj)R2 11^2 + II '^'hi,± 112,2) • 

The second part, namely the f^A-moment of the 7rr-term, is estimated as before using 
In order to control R2v^a{0 7^ we take the inner product of (12.291) with v^A , 

(-l)«^r,p +^^^y^A,R2{Q)+eG j ^{v,A)MR2{Qdv 

= ^ivzA,L7R2i^.))+g,^,^AiQ + iv,A,lUR2)iQ). 

Hence 

(-1)^" ieG f d - - i - 

R2.iA{iz) = ^^^r.,2A + / j—{vA)MR2{Qdv - —{v^A, LjR2{Q) 

-^9i,.Mz) - UvzA,IhiR2)iQ). 

Sz Sz 
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Combining this with fl2.33p and noticing the pairwise cancellation of the r-terms with 
positive and negative ^'s, gives 



\R 



■2v'}A 



< C(- II Z/5(/ - Pj)R^ \\l + II u 2^1,^ ||2 +7] II PjR2 Wlo) 



and so using fl2:27j) and fICTD 

|44(0)|' < C{\ II U-^I-Pj)R2 III2 + II U-kl,± Wy+VW PjR2 ||2,2 



< II / 



- I|2 
2,2 



"S-LX 112,2) +V II ^J^2 ||2,2 



The Vx- and Uy- moments are analogous to each others, so we only discuss the former. This 
can be treated similarly to the previous -R24(0) case but with A = 27ri?2(=F7r, ± 7rr(w), 
Vz ^ 0. The 27ri?2-term is now ingoing, and its inner product with gives zero. The 
vrr-term is estimated as before. For the sum in fl2.33p of the other Fourier coefficients, we 
notice that R2vAC) = R2vlvSC) — R2LvlvSi)i we can proceed as before. Since PPjR2 
differs from PR2 by terms of order e which are already under control, we can summarize 
the results so far as 

J (|^2^,(0)|2 + |i?2^,(0)|2 + 1/22^3(0)1^ + 14^,(0)1^) dz (2.36) 

< C{\ II Il2,2 + II ^'hl,± 112,2) +^ II R2 ||2,2 • 



Finally for the ^2o(0)-moment, start by considering (12.281) with the new boundary 
conditions 

R2i-^,v) = ri-^,v), Vz>0, (2.37) 
^2(^,^) = ^^rT / [R2{7r,w) - f-{-n,w)]wzM{w)dw, Vz < 0. 

The new boundary values are constructed in such a way that ^^^R2{—'n,v)vzMdv = 
and this property will allow the new boundary conditions to be equivalent to the old ones. 
In fact, since /^^^^ M+(t>)t>^(it> = —1, 

v^R2{7r,v)M{v)dv = / VzR2{Tr,v)M{v)dv + / VzR2{Tr,v)M{v)dv 

Jvz<0 Jvz>0 

Wz[R2{7^,w) - f'{-7r,w)]M{w)dw + / VzR2{n,v)M{v)dv 
Wzf^ (— vr, v)M{v)dv. 

«>z>0 
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Existence and uniqueness for the new problem are well known. We shall verify that the 
new problem also satisfies the old boundary conditions of (12.281) . Kt z = — vr, for f z > 
the new ingoing boundary condition for R2 is R2{—7!',v) = f^{—7r,v)), which coincides 
with the second equation in f l2.28p since v^27r ^QR2{—7i,w)\wz\Mdw = 0. At z = n the 
new ingoing boundary condition is 

R2{n,v) = M~^M+ [R2{7i,w) - f~{~7i,w)]w,M{w)dw, v, < . 

This coincides with the old boundary condition at vr, given by the second equation in 
(E2H]), that is with 



R2{7i,v) = r{n,v) + M'^M+ R2{7i,w)M{w)wJw, v,<0, 



provided that 



-f-{n,v) = M-^M+ f-{-7i,w)w,M{w)dw, v,<0, 



or, recalling the definition of / , 

[ {Ri{7i,w) + ^i^^^^)M{w)wJw = I {Ri{-Ti,w) + 'i^^^^^)M{w)w,dw. 

To check this we note that from the assumption that the inhomogeneous term g in the 
equation for Ri is such that f^^ dvM{v)g{y) = 0, it follows that 

{v^,Ri{n,v)) = {v^,Ri{-n,v)). 

Then, using the boundary conditions for Ri this becomes 

M{v)Ri{7r,v)vJv - [ M{v)^{-K,v)vjv 



+ 1 M{v)-{7r,v){-n,v)vjv- I M{v)Ri{-n,v)v^dv = 0. 



Dz>0 ^ JV2<0 

Using {ip{±TT, v),Vz) = 0, the claimed equivalence follows. Thus R2 with the new boundary 
conditions equals the unique solution to fl2.28p . 

We write f l2.28p with R2 = R2k{z). The left- hand side becomes 

-It \ ( ^-^2 n^^2 

K iz) \ v,—^ - eG ^ 



^ dz dvz , 

Multiply the equation for R'2 by MvzK,{z) and integrate over [— tt, z] x M^. The l.h.s. gives 

i^i(...)M.j*.-/i^i(-..,.)M«;*. 

-eG j dz' j R'^iz' ,v)Mvldv + eG j dz' j R'2iz',v)Mdv . 
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Since J dvvlR2Mdv = R20 + + ^ rfff^-R^Mrff , by integrating the equation again 

over [— 7r,7r] we control |-R2o(0)| in terms of known quantities plus terms bounded in the 
II ■ ||2,2-norm, multiplied by a factor and the integral J R2{—7^,v)Mv^dv. The contri- 
bution due to the incoming part is given by /~. Therefore, we need an estimate of the 

outgoing boundary term / i?2(— tt, v)Mv1dv. To do this, we repeat the steps above but 
this time integrate over {[— tt, tt] x M^; Vz > 0}, to get on the l.h.s. 



R2{7T,v)Mv;dv - / R2{-TT,v)Mvidv 

/TT /• /*7r /» 

dz R!^{z,v)Mvldv + eG j dzl R2iz',v)Mdv. 

Since R'2 solves the problem with modified boundary conditions, the incoming part in 
— TT can be bounded by ||/~||2,2,~- Therefore we get an estimate of / R2{7i,v)Mvldv in 

Jv^>0 

terms of the norm of /~ and again quantities bounded in the || ■ ||2,2-norm multiplied by a 
factor e. To obtain the estimate of ^QR2{—7i,v)Mvldv, we integrate over [— vr, vr] x M'^. 
The final estimate is 

|^20(0)| < C{- II ||2~ +- II - Pj)i?2 ||2,2 +V II R2 ||2,2 + || J^'hl l|2,2)- 

e e 

By combining with the other moment estimates (12.361) . and using (12.271) . we obtain that 



PjR2 \\l2< - II r IIL + II ^~-'9i Wh ■ (2.38) 



This completes the proof of the lemma. □ 

Based on this 1-d analysis, it follows in the 2-d case that 
Lemma 2.7. The solution R2 to /12.28\) satisfies 

II PjR2 ||2,2< c\ II /- IIL • 

Proof of Lemma 12.71 Consider equation (12.281) for the Fourier transform in x,z of R2, 

R2 = ^ 2-^2, 

i{lii.v, + izVz)R2 = -L^2 + eGM-^^!^^ + I^ilh) - v^rii,, v){-lf% (2.39) 
e dvz 

ri^x, v) = TxR2{ix, vr, v) - J^xR2{^x, -vr, v). (2.40) 

In the case 7^ 0, 7^ we can reach, as in the case of Lemma 12.61 a bound 
like (I2.3ip . If we consider the Fourier components with ^ large, we see that in (I2.3ip 



23 



the r-terms are multiplied by a small number for ^ large and then can be estimated, by 
using (12.261) . by t] || R2 \\2.2j with a small 77, plus the earlier terms. We notice that in 
the case C,x = 0, equation (12.39^ reduces to (12.29^ . Hence, we can apply Lemma [2^ to 
i?2(0, ^z, = J R2{x, z, v)dx and, taking into account that gi is of order 5, get a bound for 
the Fourier components PjR2{0,^z), for S small. The remaining components in the case 

7^ 0, ^x,(,z bounded, can be estimated in the following way. 
We start from the moment R2v:,- Notice that r„^(^x) = fv^^x), where /~ is the function 
defined as 

j-^^ ^ U^xf'){^x,'^,v) for Vz < 0, 
Hence, by integrating (12.391) we obtain 

\^x\\R2.A^x,o)\<c\i-i^x)\. 

Then, we look for a bound for R2vA^x, ^2) when 7^ 0. We use the function B introduced 
in (I2.35p . B being the solution of the equation L{vxVzB) = VxVz- 

We have {vlvxB,R2) = {vlvxB,{I - P)R2) + {4^x,Vxv'^^B){iIjx, R2), where ip j j = 
0, . . . , 4 are as usual the vectors of the orthonormal basis in Kern(L). We multiply (12.391) . 
written for = 0, by v'^B — {Tpx.VxvlB) and integrate over velocities. We first use the 
relation so obtained for = and obtain, 

KiB{ix)\ < c(i II u--2Ljif^j)R2{U0) h +eG \\ mx,0) ^ +\f-i^x)\ 

+ II U-'^lUR2m,0) h + II {I - P)R2{^x,0) II2), 

having used the fact that |r^,^| < c|/^|. We use again that relation for 7^ 0. Since 
{v^B,R2) = {v^B,{I — P)R2) + {R2,4'z){i'z,v^B), we obtain in this way an expression 
for ^zR2v^{^x,C,z), for ^2 7^ 0, in terms of quantities under control, since with the previous 
subtraction we have removed from the equation the term {ipx, -^2)- As a result, for ^z 7^ 0, 

le.l I R2vMx,Q I < II U-"2LjiI^j)R2{^x,Q II 

+ II U-^lUR2m.^z) II +e II R2{^x,Q II 

+1 II iy-'^Lj{r^j)R2{^x,0) II +e II R2{^x,0) II 

+ 11 u~-^ihiR2m,o) ii+ie.i wri^x) 11). 

Notice that the last term is bounded because ^x is bounded in this part of the proof, the 
terms with ^ large having been estimated before. 

In the same way the f^j-moment can be controlled when ^2 7^ 0, 

1^x1 I R2vA^x,Q \< C{- II U"-2Lj{f^j)R2{^x,Q II + II U"-2MR2){^x,Q II 

e 

+e II II +1 II i^--^Lj{r^j)R2iU0) II2 +e \\ M^x^O) \\ 

+ 11 u'iiuR2)i^x,o) II +i/-(e.)i). 
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For the f^-moment R2vX^xi 0), by using the proof of Lemma [2 .Gt 

\R2.Ai^M < riQ II + II J'.R2vAQ ID- 

To proceed, we observe that 

II ^x.i?2..(ex) II < C[^{J2 II ^-^^j(/^j)i?2(ex.,e.) f )^ 

+(5^ II u--^lUR2){C.,Q ft^+e{J2 II ^2fe,e. 



+- II u-"^Lj{I^j)R,{^,,0) II +£ II M^ccO) II +|/-(e. 



+ 11 u--2iMR2){^,,0) II ). 



And so, 



|i?2..fe,0)| < C( II r(e.) II +e II i?2 ||2,2 +^($^ II iy-^^Lj{I^j)R2{U^z 
+ (5^ II u-"^H^,){^^,Q f )l + i II u-"2Lj{I^j)R2{U0) II 
+ II z/-^i^(^)(ex.,0) 



To estimate R2vy{^x, 0), we use the method of Lemma [2^ that is the proof of (12.361) . Then, 
to get the estimate for 7^ 0, we multiply (12.391) by VyVz and use v^VyV^ G (Kern(L))-'- to 
estimate the term in r as 



|r.2,^(e.)| < C{- II u-Hj{f^j)R2{^.,0) II2 +6 II i?2(ex.,0) II2 
+ II u-"^H^2m,0) II2+ II (/-Pj)4(ex,0) II2). 
So with C depending on ^a,, 

II i?2.,(ex.,0) II2 < Ci-iYl II ^"^A/a'^j)^2(ex,e.) 112)^ 



e 



+{Y1 II 112)^ + II ^2(e.,e 

+ ^ II Z/-5Lj(/^j)i?2(ex,0) II2 II R2i^.,0) II2 

+ II z/-^i^(^)(e.,0) II2 + II (/-Pj)f?2(ex,0) II2), 



zj 112) 
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and similarly for R2vy{^x,^z) when 7^ 0. 
For the '?/'4-moment we shall use 

^2.^^(0 = ^^24(0 j vyAMdv + R2±,2a{0- (2.41) 
Here in the 2-d case, 

— {v,A)MR2{0dv 

The additional term in comparison with the 1-d case, belongs to R2±, and can be estimated 
by (12.271) . An estimate of the boundary term r^2A can be obtained by multiplying f l2.39p 
with VzA for = and integrating. For ^2 7^ this gives 

l^2.iA(e.,6)r < II u'^I-Pj)R2 +s' II R2 ) 

with C depending on ^. Using f l2.4ip the same estimate holds for R24.{C,x,^z)- For 7^ 
the Fourier component R2v^a{^x,0) can be expressed by (I2.33P and (12.390 . including the 
pairwise cancellation of the r(^^) terms. Treating the A-term as in the 1-d case, gives 

\R2vu(^x,o)\' < II uHi-Pj)R2 Wh +e' II R2 Wh + II riQ IIL ) 

with C depending on C,x- Again by (12.411) the same estimate holds for R2'i{^x, 0). 

The ipo moments when 7^ may now be obtained by multiplying (I2.39P with Vz and 
integrating. Arguing as above and using the earlier estimate for the ?/^4-moment we get 

\R2o{^x,Q\' < II uHl-Pj)R2 +e' II R2 + II f-{^x) IIL ) 

with C depending on C,- And so there only remains the tpo moment for R2{^x,^) when 
7^ 0. Multiply ( 12.39P for ^ = {^x, 0) with Vx and integrate. For the boundary term 
Ty^vA^x), multiply (I2.39P for ^ = (^^, 1) with Vx and integrate. All terms in the upcoming 
expression for r^^jjA^x) are then under control. This gives 

|i?20fe,0)P < C{\ II v"^{I-Pj)R2 II R2 \\l2 + II r{ix) IIL)- 

Combining all the above estimates gives the statement of the lemma, 

II PjR2 Wia— ~2 II f Il2,2.~ ■ n 

The step from to for R2 follows as in the i?i-case. These estimates together give 
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Lemma 2.8. A solution to the R2-problem satisfies 

II - Pj)R, < c(s II - Pj)g +i II Pjg +^ || ^ W^,^ ), 

II PjR2 III2 < II iy-Hl-Pj)g +1 II Pjg +1 II ^ III2,. ), 
II ^^^2 ||L,2 < II - 11^,2 +^ II III2 II ^"^^7 llL,2 

+ ^ II ^ lit)- 

The previous estimates can be used to prove 
Theorem 2.9. There exists a solution R in L|^([— tt, tt]^ x M^) to the rest term problem 

v^-VR- eGM-^^!^^ = -LR + J(R, R) + HiR) + ea, (2.42) 

ovz e 

f 1 - 1 - 

R{x,^'7i,v) = / {R{x,^TT,w) -\ — ■ilj{x,^n,w))\wz\M^dw ilj{x,^Tr,v), Vz^O, 

Proof of Theorem 12.91 The rest term R will be obtained as the limit of the approximating 
sequence {-R"}, where = and 

■ Vi?"+' - eGM-^^^^^^—^ = -LjR''+^ + H^iR^^^) + J(/?", i?") + ea, 



R-^\x,Tn,v) = ^ 



I lb lb 



Here (/ — Pj)g = e{I — Pj)a is of order four, and Pjg = ePja of order five. In particular, 
the function R^ is solution to 



■ VR' - eGM-'^^^^ = -LjR' + Hi{R') + ea, 



dv 



R^{x,TT^,v) = ^ [ {R^{x,TT^,w) + —{x,TT^,w))\wz\M^dw--{x,TTr,v),Vz'^0. 

Splitting R^ into two parts Ri and R2, solutions of fl2.23p and (12.241) with g = ea in (12.231) . 
then using the corresponding a priori estimates. Lemma [231 and Lemma [?!8l together with 
the exponential decrease of i/j, we obtain, for some constant Ci, 

II ^^R^ l|oo,2< CiE^, II U^R^ ||2,2< CiE^. 

By induction for e sufficiently small, 

II z^^-R^ ||oo,2< 2ci£:^ j<n + l, 

II Z/^(i?"+l - R^) ||2,2< 028^ II Z/^(i?" - i?"-l) 112^2, n > 1, 
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for some constant C2. Namely, 

e dv. 



Here, = (/ - P)G'"+^ = J(i?"+i + i?", i^'^+i - i?"). It follows that 



II z/5(i?"+2 - ||2_2 < ce-^ II zy-^G^+i ||2^2 

<C£-^( II Z/^i?"+l ||oo,2+ II //^i?" 1100,2) II 

Consequently, 

II ||2,2<|| I^HR"^^ - ^"^') l|2,2 +...+ II I^HR^ - R^) h,2 + II ||2,2< 2Cie'^ , 

for e small enough. Similarly || i?""*"^ l|oo,2< 2ci£§. In particular {R^} is a Cauchy sequence 
in L|^([— TT, vr]^ x R^). The existence of a solution R to fl2.42p follows. □ 

From here Theorem 12.11 follows, and as a consequence the first part of Theorem II. 1[ 
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3 Stability: the expansion. 

In the previous section we have constructed a stationary solution Fg of the Boltzmann 
equation close to the clockwise roll hydrodynamic solution hg. In the next two sections 
we study the behavior in time of a small perturbation of by writing the perturbation 
as a truncated ^-expansion and in particular in this section we show the decay to zero 
in time of the first terms of the expansion. This result relies crucially on the hydrody- 
namical stability under small perturbations of the hydrodynamic roll solution hg. Hence, 
before starting the construction of the Boltzmann solution, let us recall some known hy- 
drodynamic results. The Oberbeck-Boussinesq (0-B) equations with periodic and rigid 
boundary conditions (see [Jo]), describing the hydrodynamic behavior of the fluid in the 
present setup in dimensionless form, are: 

dtU + u-Vu = fiAu~Vp-e,Ge, (3.1) 

^{dtd + u-ve + Xu,) = Hao, 

div u = 0. 
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where is the unit vector in he positive 2;-direction, m G and G R are the velocity 
field and the deviation from the linear temperature profile respectively, 17 and k are the 
dimensionless kinematic viscosity and conductivity respectively. The initial conditions are 

u{x, z,0) = uo{x, z), div mq = 0, 9{x, z,0) = 9o{x, z) 

for any x, z & = (— /uvr, fin) x (— vr, tt). The boundary conditions for this problem are 

u{x, — TT, t) = u{x, TT, t) = 9{x, —71, t) = 9{x, TT, t) = 0, X G [— TT, Tt], t > 0. (3.2) 

Here the notations are as in the Introduction. For a proof of the existence of a global in time 
solution for small initial data see for example [FMT] . The laminar solution is the trivial 
stationary solution u = 9 = 0. It is the unique solution for Ra < Ra^. and is asymptotically 
stable for Ra < Ra^ [lu], |FMTj . After Ra^ a pair of new stationary solutions appear. In 
[lu] and [lul] it is proved that there exists 5q such that for any < 5 < there are two 
stationary roll solutions {uf , 9f) corresponding to the Rayleigh number Ra = Ra^l + S), 
of the form 

u^{x,z) = tS Co (j) + 0{6'^) (3.3) 

9fix,z) = tS Cot + 0(5''). 

Here Co is a positive constant, the couple (0, r) is the eigenfunction corresponding to the 
smallest eigenvalue do of the linearized problem, namely the solution of 

17 A0 — Vp = CzCt, div = 0, kAr = do(pz, 

(0, r)(x, -vr) = (0, t){x, tt) = 0, (0, r)(x, z, t) = (0, r)(x + /ivr, z, t). 

In |Iu2j both solutions are proved to be stable for small perturbations (see also [MWj ) . 
All the previous results are stated in the Sobolev spaces if2, but, by general theorems on 
PDE of parabolic type [La] (or by the method in [GhJ), the regularity can be improved to 
higher Sobolev spaces Hk. Hence we can state the stability theorem in a form suited to 
our purposes. Let {us,9s) G {H^Y, k large enough, be the clockwise solution of fl3.3p . 

Theorem 3.1. Let {u,9) be the periodic solution of the following equation 
dfU + Us ■ Vm + u ■ Vus + u ■ Vu = fjAu — Vp — ezG9 
^{dt9 + Us-V9 + u-V9s + Xuz) = hA9 
div u = 0, 

u{x, z, 0) = ^0(3;, z), 9{x, z, 0) = 9o{x, z) (x, z) G [— /uvr, jin] x [— vr, vr] 
u{x, — vr, t) = u{x, vr, t) = 9{x, — vr, t) = 9{x, vr, t) = 0, x G [— vr, vr], t > 0. 

If{uo,9o) G {HkY , k sufficiently large, and \\ uq \\h^ + || 9o ||_ffj.< no, for no small enough, 
then, {u,9){x, z,t) is in {Hk)^ and lim{u,9){t,x, z) = exponentially in time in {Hk')^, 

for any k' < k. 
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Notice that the convective solution hconv = {us, Tg) in the Introduction, fll.Sp , is related 
to Us, 6s in (13.31) through the shift Tg = 6s + \{z + tt). 

In the previous section, we have constructed a positive stationary solution of the Boltz- 
mann equation (11.11) as Fs = M(l + $^). Here, we want to study the evolution of positive 
perturbations of Fs. The perturbation defined through F = M(l + $^ + $'^), with F 
a solution to (11.11) . has to solve the initial boundary value problem 



9$^ 



1 G d(M^^) 1 / \ 

+ -v^ . V*^ - = ^ [L^' + ^(*^ + , (3.4) 

$^(0, X, z, v) = Cq{x, z, v), {x, z) G (-tt, tt)^, V G M^, 

$^(t, X, ±7r, v) = I \Wz\M^%t, X, ±7r, w)dw, ^ 0, t > 0, x e [-n, tt]. 

The following initial perturbations F{0,x, z,v) — Fs := $^(0,x,2,f) = (o{x,z,v) are con- 
sidered, 

5 

Co(x, z,v) = J2 ^"$^"^0, X, z, v) + e'p,, (3.5) 

n=l 

where 

F{0, ■,-,■)> 0, II P5 ||oo.2:= sup / sup pl{x,z,v)Mdv ] < c, (3.6) 

e>0 yj (x,^)e[-7r,7r]2 J 

for some constant c. The nonhydrodynamic part of the functions $'•"^(0, x, z, f ) is deter- 
mined by the expansion as explained below together with some terms of the hydrodynamic 
part. We will denote by l'f'\t, x, z) the coefficients of the functions ipi in the hydrodynamic 
part of ^("^(t, x,z,v). The functions lf'\t,x,z) will be determined by the solution {u,6) 
in Theorem 13. 1[ Finally, we require 



X 



C,q{x, z,v)MipQ{v)dxdzdv = 0. 

7r,7r]2xR3 



Since the Boltzmann equation conserves the total mass, it follows that $^ will satisfy 

M$^(a;, z, v)ilJo{v)dvdxdz = 0, t > 0. 



'[-7r,7r]2xM3 

We write an ^-expansion for $^ in the form 

5 



^''{t,x,z,v) = ^^^''\t,x,z,v)e'' + eR{t,x,z,v). 



n=l 



For the proof of stability we need to show that $'^"^(t, x, z, v) converge to zero, when time 
tends to infinity in a suitable norm. To this end we will construct explicitly the first terms 
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of the expansions. The behavior of the higher order terms will then be evident from this 
analysis. This construction is by now standard and contained in many papers. We give 
here a sketch of the argument for sake of completeness, and follow closely the analysis in 
[AEMNj . 

In the following we use the notation {h) = J^^dvh^v). By plugging the expansion 
into (13.41) . as a first condition, has to be a combination of the collision invariants 
V'i,? = 0, . . . ,4 

so that pi = ul = = \ ul = iP, = ^/f \ We require that u\ 9^ satisfy 

the initial and boundary conditions in Theorem 13.11 and, as a consequence, do not need 
boundary layer correction to the first order in e. Indeed, m. z = — vr the solution is already 
of the right type. On the other hand, M + e^^^\ when evaluated for 2; = vr, cannot satisfy 
the boundary conditions, but differs from it by terms of order e^, which will appear in the 
corrections of higher order. Hence, for n > 1 the higher order corrections are decomposed 
into a bulk term i?*^") and two boundary layer terms 

To determine the functions p^, and 0^ which give $'^i^(= -B^^-*), we consider the 
equation obtained by equating the terms of next order. Note from the previous sections 
that the stationary solution can also be expanded in e, and denote by <l>i'^'' the terms of 
this expansion. The equation which we get at next order, by ignoring boundary layer 
corrections, is 

^^^$(1) + = + J($«, $«) + J($«, $«) (3.7) 

ox oz 

It can be seen as an equation in B^'^\ whose solvability conditions give the usual incom- 
pressibility condition and the Boussinesq condition 

divM = 0, V(^i + pi) = 0. (3.8) 

The Boussinesq condition fixes p^ = —9^, up to a constant. To determine 9^ and we 
look at the solvability condition at next order in e. Indeed, once ( 13. Sp is satisfied, we can 
deduce from (13.71) the following expression for B^'^\ where denotes the inverse of the 
restriction of L to the orthogonal of its null space, 

5(2) = . v$« - J($«, $«) - J($«, $«)] + ^ jf (t, z) . (3.9) 

i=0 

The coefficients are undetermined at this point and will be partly fixed by the solvability 
condition for the equation at next order in e and the rest of them in some later step, 

^5(1) + y . vi?(2) _ }_G^{MB^^^) = LB'^'^ + J(<1>(2), $«) + J($(2), $(1)) + j($(i), $(2)). 

(3.10) 
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The solvability conditions for this equation, 

+ V . VS(^) - 1:G A(M5«)]) = 0, 2 = 0, ■ ■ ■ , 4, (3.11) 

produce the equations for and 9^. Let us fix z = 1,2,3 in (13. lip . Then the first term 
gives the time derivative of u^. The third one reduces to for z = 1,2, and to —Gp^ for 
i = 3 after integrating by parts. The term {v ® vB^"^"*) gives rise to dissipative transport 
terms and a term which can be interpreted as the second order correction to the pressure 
P2- The term J($^^-', in ( 13. 9p produces the linear transport terms, depending on the 
stationary flow. The result is 

+ ■ Vu^ + ■ Wu] + u\ ■ Wu^ = fjAul - VP2 + e,Gp^ . 

ot 

Using the Boussinesq condition we replace the term Gp^ by —G9^ + const. The constant 
can be absorbed in the pressure term that we rename p. 

Remark. There are constants (one coming from the Boussinesq condition, another 
from the pressure condition) at any order which will be determined in the end by the total 
mass condition. Since we are asking that the total mass of the perturbation is zero we can 
put to zero all the constants. 

To get the equation for the temperature, one has to look at (13.111) for i = 4. It is 
actually more convenient to replace ip^ with the equivalent ifj^ = |(f^ — 5). We have 

(^4, /i) = le\ G(^4, ^/i) = -ulG, (^^4, B^'^) = -hve' + ^u'e' + uy + ey . 



Putting all the terms together, we get 



dt 



Gu\ = -kA9^ 
' 2 



This equation has to be solved with boundary conditions ^^(±l,t) = for t > 0, and an 
initial condition which is completely arbitrary. 

Remark. In the previous equation there is a term —Gu\ which does not appear in 
the usual 0-B equations. This term can be absorbed by changing the boundary con- 
ditions. Here 0\.,u\ are the hydrodynamic terms of first order in e in the expansion of 
and hence they coincide with Us,Ts in (11. 3p . The boundary conditions for Tg are: 
Ts{x, —7i,t) = 0,Ts{x,n,t) = 2\n. The shifted temperature Tg = — Gz will satisfy 
the usual Boussinesq equation, in which the term Guz is missing and a different boundary 
condition, Ts{x,7T,t) = (2A — G)n. This aspect was discussed in |EML] . It was pointed 
out that starting from the compressible Navier-Stokes equation or the Boltzmann equation 
in the scaling we are considering, one obtains a set of equations which differ from the 
usual 0-B ones for this shift in the boundary condition for the temperature. By scaling 
the variables, this amounts to the usual 0-B equations in dimensionless form, with a new 
Rayleigh number given by Ra{l — G). We conclude that nothing changes in our analysis. 
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To summarize what we got so far, {u^, 9^) has to satisfy the 0-B equations in Theorem 
13.11 By fixing the initial conditions so that the assumptions of the Theorem are satisfied, we 
get that (m^, 9^) vanishes exponentially in time, with its spatial derivatives. Since 9^ differs 
from by a constant, which can be taken as zero, we may conclude that for q G [1, +cxd], 
II ^'^^^ ||(j,2 is finite and converges to zero exponentially in time. 

The second order term in the expansion, <l>'^^\ is not yet completely determined. Equa- 

tion (13.111) with i = gives = div fixing div /*^^\ Moreover, a combination of 

(2) (2) 

Iq and II contributes to the pressure p which is determined by the previous equations, 
so that these parameters are not independent. 

The nonhydrodynamic part of S'-^^ is a linear function of the derivatives of p^, 9^ which 
are in general different from zero at the boundaries. Therefore the non hydrodynamical 
part of -B*-^-' is completely fixed (even at time zero) and violates the boundary conditions. 
We need to introduce to restore the boundary conditions by compensating the non 
hydrodynamical part of 5^^^ which is not Maxwellian. We explain how to find the correction 
b_ . The correction b^'^ is found in a similar way. Here = 2Ar^Q{M_, M_ ■). We 

(2) 

choose b_ by solving, for any t > 0, the Milne problem for z~ > 0, 



Vz^h-e^—G~ — {Mh) = L-h, {vzh):= dw^h = 0, (3.12) 
dz- M dvz Jm.3 



where z~ = e^^{z + vr) is defined as the rescaled z variable near the bottom plate, and 
G~ is a smooth force rapidly decaying to zero far from the bottom plate. Indeed, the 
gravity force has been decomposed in three parts, a force constant in the bulk and two 
boundary parts (see [EML] . |ELM2] for details). We impose the boundary condition 
at 2~ = in such a way that the incoming flux oi h dX z = — vr, Vz > 0, is given by 
{I—P)B^'^\—1, V] t). The results in |CME] tell us that as z~ +oo the solution approaches 
a function q_ {v, t) in Kern L~ . Note that in q_ there is no term proportional to ips because 
of the vanishing mass flux condition in the direction of the z axis (vzh) = 0. Thus we set 
b^^\x, z_, V, t) = h{x, z_, V, t) — q^^\x, v, t), which will go to zero at infinity exponentially in 
z~ . This produces a term b^^\x, 2TTe~^, v, t) = ip2,e{x-, vr, f , t), exponentially small in on 
the opposite boundary. Scaling again to the variable z, the resulting term in the expansion 
is thus = i?'-^-' + b^^ + 6_ , and is such that in z = — vr, for example, it has zero non 
hydrodynamic part, while the hydrodynamic part is 

4 

^^^\x,-TX,v-t) = ^lf\x,-Ti-t)i)i{v) + bf{x,2e~\v,t)-q'^':\ i;, > 0, t > 0. 

We are not yet done since M$*^^''(x, — vr,f) is not Maxwellian for Vz > 0, (as it should, in 
order to satisfy the boundary conditions) because of the presence of terms proportional 
to t/ji, i = 1,2,4 in q^^^ and b^^\t,x,2e^^). The latter is not important because it is a 
correction exponentially small in e and will be compensated in the remainder. The former 
will be compensated by the coefficients , i ^ 0, 3, that can be chosen arbitrarily on the 
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boundaries. Moreover, we have to choose /g^^ = on the boundaries, because {vzq^-) — 0. 
We are left with 



$(')(i,x,±7r,T;, ^ 0) = ajM± + ^2,.(a;,±7r,^,t), at = I^i^^) " (9± ), 

where ?/'2,e are terms exponentially small in e. Finally, we impose the impermeability 
condition (v^^^^^) = by choosing 



^ M 



I VzM[^^'^\t,X,±Tl,v) -ij2,e{t,X,±TX,v)]dv, ^ 0, t > 0. 



The coefficients = 1,2,4 of the hydrodynamical part of B^"^^ are determined by 

the compatibility condition for the equation at next order in £, 



where 

Idt M dv,^ ' ^ ' 



■'z 

4 



-J($«,S(^))-J($«,Sf))J+5]^.7f, 

1=0 

together with the boundary conditions II = {Q- )%, i = 1,2,4. Then /q is found up to 
a constant that is chosen so that the total mass associated to vanishes. Proceeding 
as in the determination of the Boussinesq equation, we find now a set of three linear 
time-dependent nonhomogeneous Stokes equations for 

dtp^ + div + div {p^v}) + div {u^) = Nq 

dtu^ = ■ Vk^ + ■ Vu^ = fjAu^ - + Gp^ + Vdiv + N 

dt9^ + ^[div + {p" + ^')div u^] + p^dtO^ + ^div u"] = ^^[A^' + {Vuif] + N^, 



where A^o, -^4 depend on the third order spatial derivatives of p\ 0^ and depend on the 
third order spatial derivatives of v}. We remember that is determined by p which 
has been found at the previous step by solving the 0-B for v},6^. On the other hand, 
= p'^ + 9'^ + p^9^ allows us to eliminate p^ from the previous equations. Replacing 
div (m^) as given from the first equation in the last one, and using the condition on P^, we 
get a set of two coupled equations for 6*^,^^ with a constraint on divw^. P^ plays the role 
of Lagrangian multiplier for this constraint. The nonhomogeneous term is controlled by 
the results at the previous step and hence is known to decay to zero exponentially in time 
in the right norms. Then, general theorems for the Stokes equation assures the existence 
of a solution for the chosen boundary conditions, vanishing exponentially in time. 
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Once B^"^^ is completely determined, the last equation gives the non-hydrodynamical 
part of 5*^^^ As before, we introduce the terms b^^ to compensate (J — P)B^^^ on the 
boundaries z = ±7r. The term is found as a solution of a Milne problem with a source 
term, which depends on the previous boundary corrections b^2^ and The procedure 

can be continued to any order. 

We notice that (J — P)$(") at time zero are not arbitrary, since they depend on <|>("^^) 
and its derivatives. We can instead assign at time zero = 1,2,4. Notice that the 

rest term R at time zero is of order e^. By using the results in [CME] and the exponential 
decay in time of we can state the following theorem 

Theorem 3.2. Assume that at time zero, for some suitably large k, 

II Ma'=/f^(0,a;,2) ||l2< oo, i = 1,2,4, n = 1, . . . 5 , 

where denotes any space derivative of order k. Then, it is possible to determine the 
functions n = 2, ... ,5 in the asymptotic expansion satisfying the boundary conditions 

$(")(t,x,T7r,t;) = [ |w,|M[$(")(t,a;,T7r,u;) -^„,,(t,x,T7r,«;)]ci«; 



the normalization condition / M^^'^^dvdzdx = 0, t G R^, and 

II '^'^"^ ||2,2,2< OO, II ||oo,oo,2< OO . 

Here, 

II / ||2,2,2= ( / / / \f{s,X,Z,v)\''M{v)dsdxdzdvy, 

II / ||oo,oo,2= sup ( / sup \f{t,x,z,v)fM[v)dv] . 

t>o ^ Jr3 (x,z)£nn ^ 



4 Stability: the remainder 

We remember that the solution to (13.41) . is written as $^ = X]i=i + 

previous section we have constructed the terms and shown that they decay to zero in 

suitable norms. In this section we construct the rest term i?, solution of 

i?(0, X, z, v) = Ro{x, z, v) = e'^p^^x, z, v), 

X, =Fvr, t;) = ^ f x, =Fvr, w) + -(t, x, =Fvr, w))|w;2|M(iw - -(t, x, =Fvr, f ), 

X G [— vr, tt], t > 0, Vz > 0, 
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Here x, ivr, t>) = x, ivr, t>) is the Knudsen part of the asymptotic expan- 

sion from {t,x,^7T,v), exponentially small when evaluated at (t, x, ivr, f ). A contains all 
the terms fully coming from the asymptotic expansion, 

where $ = Yli ^^■^''s-^ , and J p^Mdxdzdv = 0. We shall require that the initial value of $^ 
is close to zero, and in the sequel introduce smallness assumptions. 
The following norms will be used, 

II -R ||2i.2,2 ~\ / / R'^{s,x,z,v)M{v)dsdxdzdvj, 

II -R II 00,2,2 = sup ( / / / R^{t, X, z,v)M{v)dxdzdvj , 

f 1 

II -R II 00,00,2 = sup ( / sup R'^(t,x,z,v)M{v)dv] , 



t /-TT /• ,1 

2 



/ ||2t,2,~ =( / / / vM{v)\ ]{s,x,--n,v)f dvdxds] + 



Vz I M{v) I /(s, X, TT, t>) p dvdxds ) ^ , 



f ||oo,2,~ = ( sup / / t;2M(^;) I f{t,x,-7i,v) \^ dxdv) + 



(sup / I 



<0 



Vz I M(t') I f(t,x,n,v) p dxdv^ 



We will prove the existence of $^ and the stability result (11.7p . We follow closely the 
approach in Section 2, starting from dual, space-periodic solutions to a linear problem 
(in the rescaled time variable r = e~^t) discussed in the following lemma. We use the 
notations introduced in Section 2, Lj = L( ■ ) + eJ{q, P ■), but here the function q has the 
expression q = £:^^(<l> + $s) which is also time dependent. 

Lemma 4.1. Let ip{f,x, z,v) be solution to 

|£ + , + - ,GM-^ = ii}^ + g, (4.2) 

or ox oz ovz e 

periodic in x of period 2tt, with zero initial and ingoing boundary values at z = — vr, tt, and 
g x-periodic of period 2n. Set (p = Lp— < ip >= Lp — (27r)~^ J ipdxdz. 



36 



Then, if e < Eq, 6 < 60, for Eq^Sq small enough, there exists rj small such that, 

II ||oo,2,2< c(e^ II Z/~5(/ _ P)g ^.^^^ +£-3 || Pg ^^^^^ J^r]E^ ||< P^p >||2,2 ), 
II Z/5(/ - P)ip ||2,2,2< c(£ II Z/"5(/ _ P)g ^^^^^ + || Pg ^.^^^ +r]e \\< Pip >||2,2 ), 
II l|2,2,2< II Z/"5(/ - P)g) ||2_2,2 || Pg ||2,2,2 +?? ||< >||2,2 )• 

Proof of Lemma 14.11 A variant of the method in [Ma Sen 7.3] can be adapted to the 
present setting with a force term, to obtain the existence of a solution to (14.21) . 
Denote by ip{f,^,v), ^ = {^x,^z) ^ the Fourier transform of ip with respect to space, 
and define g analogously. Then for ^ 7^ (0, 0), 

or e ovz 

Here v^" = (/it;^;, t;^), r = ^^(/^(f , vr, t;) for > 0, r = J^xLp{f,^x, -7i,v) for < 0, 
with JF^, denoting Fourier transform with respect to the x-variable. Let (3 he a truncation 
function belonging to C^(M) with support in (0, 00], and such that (3{f) = 1 for f > Tq for 
some To > 0. Let = 0(3. Then, for (0, 0) ^ ^ G Z\ 

OT e ovz or 

Let JF be the Fourier transform in r with Fourier variable a. We put 

Z = J^ (e-^L)^ + g(3- \v,\rP{-l)^^^ , Z' = T {e-^L)S^ + 9^), U = {ia + i^ ■ v^")'^ . 

Let X be the indicatrix function of the set {f ; \ a + ^ ■ v'^ \< a|^|}, for some positive a to 
be chosen later. Similarly to Section 2, the elements ipo,---,ip4 are an orthonormal basis 
for the kernel of L}. Let Cs{v) = {1+ \ v \y. For ^ ^ (0, 0) 



P{X^)\\ <cY, [ xH(^,^,v)ijjMdv \\^P, 



4 

< c II c-sx$ II Yl II ^^'^^ II- II ^-'^'^ II • 

i=o 

Use this estimate on the support of x for « =11 C~s^ W^^W C-sZ' \\ . 

As in Section 2, the previous estimate also holds with respect to supp xi where the 
indicatrix function xi is taken for a = y/Si. We fix 61 so that c^/5i << 1. Then the above 
estimate gives that the hydrodynamic P-part of the right-hand side, || P(xi$) ||, can be 
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absorbed by || P{xi^) \\ in the left-hand side. The estimates hold in the same way when 
Xi is suitably smoothed around For the remaining (1 — x)(l — Xi)$ = x'^Xi^ we 

shall use that $ = -UZ. Then 

II Px'xl<^ ir< c( II Csx'xiu f + II Cs+2x'x0 f ) II C-sZ' r +© 



where 

e:=-2^ / ^iX'x\U:F(^eGM-'^^^^ + ^^)Mdvn ^^^^ 

j=0 z J 

We again replace o; in the denominator by || C-s^ ||~^|| C-sZ' \\. That gives 



p$ ir< c(|| C-.$ II II C-sZ' II + 11 -^^V^f^ 'I' + ^/^ II C-.(/ - P)$ in + e. 



<^llel 

Hence, 

T./W\3r II 2 



p$ 11^ < c((|| p$ II + II C-.(/ - 11) 11 C-sZ' II + 

II c-.(/ - p)$ ir)+©- 



<^iiei^ 



Consequently, 
II P* 11^ 



< c( II (.-sZ' f II ^_^(/ _ nil 

+ II c-s{i-p)^ ir) +e>- 



We next discuss the term ©. The first term in the first integral can be bounded by e times 
an integral of a product of M, 1 + ji^^j, a polynomial in v, \ J-(p \ and JJ or C/^. So this 
factor is bounded by eC || $ ||. And so, 

II p$ f< c[ II c-.^' ir + " ^"^f II - II' ) 

4 

-2^ I i,^X'xltj{T(p^-P^Mdv(^ I i^^x''xlWP-UZ)Mdv)\ 
Therefore for ^ ^ (0,0), 
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{P^f{a,^,v)Mdvda 
<C(1 j da \\C-s{v)TLyp)ia,^r)f + j da[\\ (/ - P)<|.(a, ■) 



3=0 

Sending tq to zero implies that 



df j {P^f{f,i,v)Mdv < C(^^j^df(^ II C-s{v)L*j^){f,i,-) f (4.3) 
+ II (/ - P)^(f , e, ■) f ) + df J u-'g\f, e, v)Mdv + 



Taking 5 and e small enough and summing the previous inequality over all 7^ ^ G Z^, 
implies, by the Parseval inequality, that 



< c 



(PipY^T, X, z, v)MdvdxdzdT 
(~2^ J ^{{I - P)vf{T'jX,z,v)Mdvdxdzdf 



00 



+ j j ^ 9 iT,X,Z,v)MdvdxdzdT+ \\ 7 (p ||2oo,2,~ +V II V Il2,2,2 

As in Section 2, to use an argument based on a variant of Green's formula, we multiply 
the equation (14.21) by 2(f>MK, and integrate over [0,T] x [0,27r]^ x M^, integrate by parts 
and obtain, by using the spectral inequality (12.81) and the bounds 1 < k{z) < e^^'-^'^, 

II \\lf,2,^ + II f \\It,2,2 +1 II '^H^^ - P)V Il2f,2,2 

< C{e II I^"5(/ _ P)g ||2^^2,2 +Vl II \\lfX2 +^ II ^9 Il2f,2,2)- 

Inserting this into the previous inequality, the lemma follows. □ 

We next decompose the operator H in the remainder equation in accordance with the 
operator Lj: H{-) = J{q,P-) + Hi{-). We notice that Hi{R) is of order zero in e, 
and only depends on the nonhydro dynamic part (/ — P)R. As in Section 2, to solve the 



39 



equation for R we shall use an iteration procedure based on the decomposition of R in the 
sum Ri + R2, where Ri and R2 are solutions of two different problems. -Ri solves 



dRi 1 „ 1 dRi G d{MRi) 1 „ 1 
i?l(0, X, z, v) = Rq{x, z, v), 

Ri{t, X, =F7r, v) = -^i'ii, X, =Fvr, f ), t > 0, ^ 0, 

Here Ri is periodic in x of period 27r, and g is some given function, x-periodic of period 
2n with f Mg{-, x, z, v)dxdzdv = 0. For the existence of the solution, see the discussion of 
(14. 2p . The equation for R2 will be introduced below in (14. 7p . 

The non-hydro dynamic part of -Ri is again estimated by Green's formula; multiply (14. 4p 
by 2RiMk, integrate with respect to the variables {f,x,z,v) over [0,T] x [0,27r]^ x M^, 
integrate by parts and use the spectral inequality for Lj and the bounds 1 < k,{z) < e^^*^'^, 
to obtain, for every r^i > 0, 

II I'Rl \\lf,2,^ + II ^l(^) Il2,2 +1 II - Pj)Rl Il2f,2,2< c( II i?0 ||2,2 (4-5) 

+e II Z/"^(/ - Pj)g ||2y_2,2 +y II PjRl Il2f,2,2 +^ II \\It,2,2 +^ II ^ Il2f,2,. 

An a priori bound for PjRi is obtained in the following lemma based on dual techniques 
involving the problem (14. 2p . Consider first the problem (14. 4p without the term Hi{Ri). 



Lemma 4.2. Set h := PjRi. Then 

II h \\l2,2< C(|| Ro ||2,2 + II ^'hl - Pj)9 Il2,2,2 +^ II Pj9 Il2,2,2 +^ II Il2,2,~)- 

Proof of Lemma 14.21 In the variables (f,a;,2;,f), the function Ri is 27r-periodic in x and 
solution to 

dR. dR, OR, ^^^_, d{MR,) .... 
-K- + ■ ^ + ■ sGM — = -LjRi + g, (4.6) 

OT OX OZ OVz e 

Ri{0,x, z,v) = Ro{x,z,v), 

Ri{f, X, =Fvr, v) = —i'i'^, X, =F7r, v), r > 0, ^ 0, 

Let (y9 be a 27r-periodic function in x, solution to 

9v> dip dip .d(Mip) 1 ^ 

OT OX Oz ov. e 
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with zero initial values and ingoing boundary values at z = — tt, vr. Multiply the equation 
for <f by kMRi and the one for R by nMip sum them and integrate by parts. Then, 



dfdxdzK{z) — {Ri,ip)H + / dfdxdzdv f M—(vzkRilp) — eGk — - — 

dr J \ dz 

= j dfdxdzMKdv[^{Lj{{I - Pj)Ri){I - P)^) + \{I - Pj){Ri)L*j{I - P)^) 

+g^ + hPjRi]. 



This gives 



l2f,2,2 < -y II ' ■> ■) Il2,2 +^ II Vi'^ ) •) Il2,2 



Ki „ _ „2 1 „ _ „2 

+ Y' II ^ ^1 Il2^,2,~ +^ II 7 Il2f,2, 



+ ^ II - P^)i?i ||2 +^ II - P)^ III 

Ka .. 1 _ „ . „o 1 



+ ^ II '^'HI - P.j)9 Il2f,2,2 +^ II ^^(^ - ^)</' Il2^,2,2 

II i II 2 -^2 II 2 II 2 '^'^ — II 2 

^''V Il2,2 +— II Pj9 Il2f,2,2 +777^ II Il2f,2,2 +7r7^ II '^"'Z' Il2,2, 



for any positive constants Kj, j = 1, ...,4. All the terms computed at time r on the l.h.s 
can be estimate using Lemma [4. II and fl4.5p . leading to 



h\\l2,2< c[{K, + K,) II Ro \\l 

+ + Il2,2,~ +(^^1 + ^4 + i^3£) II ^^-^(/ - Pj)g 11^,2,2 



II 2,2 

-— + —- + -- + — + ^ /I ,2,2 + — + — + ^2 Pjg ,: 



+ II PjRi Il2.2,2 (r/ii^i +r?ii^3) + ^ + ^ + 1^ + ^) ||< Pj^if >\\U 



K, Ks K, K; 
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We are left with the term < Pj*(p >. As in Section 2, we use an approach based on ordinary 
differential equations for the Fourier transform with respect to the time and x-variables. 
Namely, the quantity < ip >:= (27r)~^ J ip{-,x,-)dx satisfies a 1-d problem including a 
small perturbation of magnitude 6 from the value at the bifurcation point. After a Fourier 
transform in time (Fourier variable a) the case of lea | < o"o with < (Tq sufficiently small, 
can be handled as in Lemma [2.41 For the remaining ex's use the term ieaTfTxP^{o', z) 
to express the ipQ-TdOTdeYit. With C2 = {v1,'ip4){vlA^ip4)^^ , project the equation along 
Vz — C2VzA, and along ipQ and use the equation, leading to an expression for 

d d d d 

^f^x{e-g^P{^v,-C2y^v,A) + -g^{Py^o + ^2)) and J^fJ^x{e—{PyDo) + —{P^vJ), 
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with '§2 a nonhydrodynamic moment of <y9, thus to an expression for 

(For details cfr Lemma 5.5 below). Estimating the ingoing terms, this results in 

||< Pip >||2,2 < C ||< Pip >x||2,2,2< C(||< Pip >x||2,2,2 || ^ ||2,2,2) 

<-\\<h >a; 2,2,2 +ri l^'"^ 2,2,2< " 2,2,2 I^^y^ 33,2 • 

So choosing £ < 1, then Ki and i^^s (resp.i^'2) of order (resp. £~^ ) and r/i of order £, 
leads to 



h ||2,2,2 < II i?0 ||2,2 + II ^ Hi - Pj)9 Il2,2,2 +^ II ^J^? Il2,2,2 



-£ 

+ ^ Il^ll2,2,~ +^ ll<^J^l>ll2,2,2 

This ends the proof of Lemma [4.21 when coming back to the t- variable. □ 

We now give the final estimates for Ri. 
Lemma 4.3. The solution Ri to satisfies 

1 „ „ „ 1 



II ^^^Rl ||2,2,2 < c( II Ro ||2,2 + II V^(/ - Pj)g ||2,2,2 +" || PjQ ||2.2,2 H 3 

V £ £2 

II -Rl ||oo,2,2 < c( II Ro ||2,2 + II V^(/ - Pj)9 ||2,2,2 +" || Pj9 ||2,2,2 H T II '</^ l|2,2,. 

V £ £2 

II l^^Rl ||oo,oo,2 < c(^£"^ II Ro ||2,2 + || Rq ||oo,2 +" || - Pj)g ||2,2,2 

1 „ „ „ _i „ _5 „ - „ 1 , 

+ ^ Pj5( 2.2,2 00,00,2 M 2,2,~ +- I W ||oo,2,- 

£^ £ 

Proof of Lemma 14.31 The solution Ri of (14.41) without i^i-term satisfies 

^ II l~Rl ||2,2~ +SUp II Ri{t) II2.2 +- II - Pj)Rl ||2,2.2< cf II i?o ||2,2 

+ 4 II ^ l|2,2,~ + II - Pj)9 ||2,2,2 +^ II PjRl ||2,2,2 +^ || PjQ ||2,2,2 

£2 ''7v^ 



for any > 0. Moreover, it follows from Lemma [4.21 that 

II PjRl ||2,2,2 < II i?o ||2,2 + II Z^"^(/ - Pj)g ||2,2,2 +^ || ^Jfi- ||2,2,2 
1 



+r77 II 112.2,' 

£V£ 
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Choosing rj = ^/e leads to the first two inequahties of Lemma 14.31 Then, to get the 
estimates, one has to study the solution along the characteristics. This analysis is 
complicated by the presence of the force, but can be done along the lines in [AEMN] and 
the result is 

II l^^Rl ||oo,oo,2< c(^^ II Ri ||oo,2,2 + II ^(0, ■ ) ||oo,2 +£ \\ 9 ||oo,oo,2 + || 7"^^! l|oo,2,~ ), 

which leads to the last inequality of Lemma 14.31 Adding the term e~^Hi{Ri) does not 
change these results. □ 

The remaining part R2 of R satisfies the equation 

dR2 dR2 dR2 i9(Mi?2) 1, „ 

i?2(0, X, z, v) = 0, 

R2{t, X, =F7r, v) = ^^jj-r- [ (Ri{t, x, =F7r, w) + R2{t, x, =F7r, w) 
H — ■ip{t,x,^7i,w)^\wz\Mdw, t>0,Vz^O, 



Its analysis is more involved and requires a careful study of the Fourier transform of R2. As 
with the stationary case in Section 2, existence for the problem (14. 7p can be adapted from 
the corresponding study in [Ma], if one includes into that approach the spectral estimate 
for Lj, and the characteristics due to the force term. 
In ( 14. 7p the given indata part is 

f~{t,x,TT^,v) = ^ [ (Ri{t,x,TT^,w) + -tjj{t,x,TTT,w))\wz\Mdw, i;^ ^ 0, 
By Green's formula for (14. 7p . and noting that H\{R2) only depends on (/ — Po)R2, we get 

e\\R2ml2+ II I-R2 \\lt,2,^ +^ II l^Hl - Pj)R2 \\lt,2,2<\\ 1^R2 11^2,^ ■ (4.8) 

This estimate is not yet final. In fact, compared with the analogous estimate for ip and 
the boundary terms here are different, due to the diffusive boundary conditions for R2. 
We follow the reasoning in the stationary case for (12.261) - (I2.27p . Taking into account the 
differences, like dependence on time, we get 



erj 

1 

The hydrodynamic estimates for R2 are obtained similarly to the stationary case. We 
again start with the 1-d (x-independent) case, with an inhomogeneous term gi which will 



^11^2112,2^ + 7 II '^Hl~Pj)R2 \\ltX2< — II r llL,2,~ +Cer] II PjR2 ||L,2,2, (4-9) 
I-R2 ||L,2,~< ^ II r 111,2,^ +C II PjR2 ||L,2,2 • (4.10) 
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take into account the x-dependence in later proofs. Reduce the equation (14. 7p to a 1- 
d problem for J dxR2{t,x, z,v) := R2{t,z,v), with an inhomogeneous term gi such that 
910 = /ffi3 dvgi = 0: 

+ - eGM-^^-^ = hjR2 + H,iR2) + g,. (4.11) 

at oz ovz e 

Then, for the solution of this new problem the following lemma holds: 
Lemma 4.4. 

II PjR2 Il2,2,2< ^ II r Il2,2~ +C2(|| PjRl ||2,2,2 + II ^'~'9l Il2,2,2)- 

Proof of Lemma 14.41 The equation for the Fourier transform, R2 of R2{t^ -2, v) with respect 
to the space variable z is 

e^^R2 + ^v,^A2 - eGM-'^{MR2) = e-^L^2 + MR2) - v-^-lf^ + 9i, (4.12) 

r{v) now denoting the difference between the ingoing and outgoing boundary values, 

r{v) = R2{t, vr, v) - R2{t, -vr, v). (4.13) 

Starting from the method of Lemma HIT] with Lj instead of its adjoint, and considering 
the ingoing boundary values as known, we reach (14.31) for 7^ with obvious changes, 

/CO n -i roo 

dt J \PjoR2\\t,^,,v)Mdv < C[^J^ dt{ II C-siv)LjR2it,^.,-) f 

dt j i-\-^\g,\\t,i,,v)Mdv 

foo II /K, 1^ II 2 



Here, similarly to Section 2, we have fixed z = Zq, t = to in the basis elements t/'i, '?/'4, 
writing Pjq for the corresponding kernel projection. At the end we replace Pjo with Pj, 
since Pj — Pjq = 0{e). The sum R = R1 + R2 satisfies by hypothesis (-R(0, ■ )) = -Ro(O) = 
for t = 0. The coefficients in the asymptotic expansion can be chosen so that, moreover, 
Ro(t, 0) = for t > 0. That gives an estimate for -R2o(0) in terms of Riq. 

For an estimate of the r-term, for z = vr we use (14.101) giving an estimate of the outflow 
at TT in terms of known quantities and PjqR2. Then, for > ^z, and large enough, the 
latter appears multiplied by a small factor and can be absorbed by the total sum on the 
left-hand side, taken over ^z ^ ^z- Hence, for ^ large the bound stated in Lemma 14.41 is 
proved. 

For the remaining bounded number of ^'s, each hydrodynamic mode is estimated sepa- 
rately. We obtain estimates for all the hydrodynamic moments in terms of ^20; which will 
be the last one to be estimated. 
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For the f^-moments multiply fl4.12p by Vz (resp. one), integrate with respect to Mdv, 
and multiply with _R2o (resp. R2vJ- Use the notation [/]_ = — 1)^^/(0). Combining 
the results, it implies that for each fixed 

We integrate over t G [0,t], to get an estimate for the term i^z[R2vAiz)\'^ on the right- 
hand side. The time derivative, then, produces a term eR2Q{iz)[R2vA*- computed at time 

t which is bounded as C£(|^2ofe)P + |^2,;.fe)P + \R2vA'^)?) < ^ll^^2||2,2(^- The latter 
is estimated by using Green's formula for (14. lip , getting 

e\\PR2\\l2 < k\\ r lib,- + II '^'hl Il2f,2,2) +^ II PjR2 \\lt,2,2 ■ 

The other terms can be easily estimated, but for the one containing the boundary r-term, 
which is equal to R2vA^z){,R2vX'^i ^) ~ R2vX^i ~^))- This term can be estimated by 

j HR2.f + -^\\ r wDdt. 

In conclusion, we have that for 7^ 

j \R2.Aiz)?dt (4.14) 

<c j dt(^\ R2oiQ p + I R2vAo) I' +v II PR2 III2 +^ II r IIL + II ^"^^71 112,2)- 

We then compute -R2„^(t,0) for C,z = 0. Multiply fl4.1ip by M, integrate over velocity and 
over z' G [— vr, 2;], followed by an integration over [— 7r,7r]. We get, after multiplication by 

R2v^ (t, 0), 



£4.,(t,0)^ + (i?2..(t,0))2 



/TV PZ f* 

dz dz'dvMgi + 2TrR2-,At,0) dvVzMR2it,v,-Tr), (4.15) 
■TT J —77 J 

where R^ := J^^dz J^^dz'dvMR2. Multiply IKT2^ by Mvz, integrate over velocity and 
multiply by R20, 

eR2o^R2v, =R2o(ry2+ J VzgiMdvJ + eGR2o j Vz-^MR2{0,v)dv. 
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Summing the last two equations, 

/TV pz n 

dz j dz'dvMgi + R2^2'K I dvv^MR2{t,v, -n) (4.16) 
-71" J —TV J 

We use this relation to bound the time integral of R^v^- We integrate over time, and use 
inequality (14.101) to control the boundary term r„2. The result is 

J \R2vAt,0)\^dt<C(^ J dt||i?2o||2 + ^ II PjR2 ||2,2+ II ^^'^l Il2,2+ II T IIL)- 

By Parseval identity, II-R20II2 = Xl^^^o l-^2o('C2) 1^ + |-R2o(0)p. The last term is equal to 
|-Rio(0)P because / dzdvR{t, z,v) = 0. The sum for > was estimated above. In 
conclusion, 



J \R2vM\'dt < C( 5^ \R20{^z)\'+ II PjRl 11^,2,2 +r] II PR2 11^2,2 + II '^hl 



|2 

I 2,2,2 



+ II - P)R2 11^2,2 +-2 II - Pj)R2 ||2,2,2 + II /" Il2,2~ )• (4-17) 



For the control of Wa;-moments, we use R2V.AC) = R2vivAC) ~ -^2j;^«2(0 (recall that 
/ vIv^Mdv = 1). Multiply (I4.12p by v^v^ (resp. v^v"^), integrate with respect to Mdv, and 
multiply with R2v^vi (resp. R2v^v^)- Adding the results implies for each that 

d ,^ . ,^ .0 ^ ^ f / A. OMR. 



dv,. 



^^{R2v^v,R2v^vl) - iiz\R2v^vl? + i^zR*2y^r3R2v:,v, j [RIv^vI^xVz 

OMR* \ f /I — \ 

+R2v,v,vy,—^jdv + Rl^^^2 j Mv^v,(-^LjR2 + Hi{R2) + v,r{-lf^ + ^ijdv 

+^2.... j Mvy,(^^LjR2 + H^2) + v,r{-l)^^ + gijdv. 

We do not estimate directly the terms involving higher moments of the boundary term r. 
To remove these terms we subtract the same expression for ^2 = multiplied by (— 1)^"^, 
getting 

£-Q^{[R2v^v,]-[R2v^v^]-) = i^z[R2v^v^]-R2v^vl{^z) + i^zR*2^^^3[R2v^v,]- 

+eG J !^[R2.o^.2Y_v^v,^{M[R2U]) + [R2,^^^^^^ 

+ [R2..v^]- J Mv,V,(^[LJR2]- + [IUR2)]- + [^l]-)^^^ 

+ [f?2...J- j Mvy,(^^[LjR2]-] + [H^2)]- + m-Ydv. 
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Integrate with respect to t to obtain for ^ 7^ 

j \R2.AQ\'dt <cj [\R2.A0)\'+v\\PR2\\l2 + \\^ki-P)R2 111,2 

+ II ^'hl Wh +^ II - Pj)R2 ||2,2 )dt- 

We have used that only (/ — P)R2 contributes to R2v^v^ and R2v^v3- 

Now we discuss the estimate of R2v^ for = 0. To ehminate the outgoing boundary 
terms, we multiply (14.111) by Mf ^ and consider first the equation we get by taking the in- 
tegral f^^ dz J^^ dz' /jg3 dv and then the one we get by taking the integral 27r dz dv. 
By taking the difference of the two equations we get for the l.h.s. (but for the force terms) 

^^tij^^j j dqVxVzMR2{t,q,v) -2^ j v^VzMR2{t,0,v)dvj + R2v^^2{t,0) 

The remaining ingoing boundary terms are zero. Now, notice that VxVzMR2{t, 0, v) = 

C4..(t,0) + /22±(t,0), where C = l^^^vlv.Mdv and 4x(t,0) = l^^^VxVzMR^{t,0,v)dv 

depends on the nonhydro dynamic part of R2- We multiply by R2y^(t, 0) and, by using the 
equation for it, we get in the l.h.s (but for the force terms) 

where T> : j J^^ J^^VxVzMR2(t,q,v)dqdzdv — 27iR2±(t,0) is nonhydrodynamic. We use 
again R2v,{0 = R2v^vl{C) - -Ri,,t,2(0 and estimate 7"i?2 by (I4.1(JI). This gives 

j \R2,M?dt <C j dt(ri II PR2 \\l2 + II l^Hl - P)R2 ||2,2 

+ 1 II z/i(/ - Pj)R2 + II /" Il2,^ + II ^-'^91 Wh ) ■ 
The f j^-moments are analogous, hence for all S,z 

j \R2.M^)?dt + j \R2.,{Q?dt <C j{'n\\PR2\\l2+\\r \\l^ + \\u'-{I-P)R2\\l2 

+ J_ II (J _ Pj)R^ + II ^-\g^ \\l^^)dt. (4.18) 



1 „ 1 

Consider next the '?/'4-moment for 7^ 0. Multiply (14.121) by v^A (resp. vlA), and 



integrate with respect to Mdv. Similarly to the proof of Lemma [2.61 this gives 

d f d 

e-{vzA,R2{^z)) +{-^)^'r,2A + iUvlA,R2{^z))+eG J —ivzAM)R2{Qdv 



^{vzA,LjR2{a) + {vzA,lUR2mz)) + {vzA,g,] 
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Sg-^ivlA + tUv'A 4^)) +eGj —{vlAM)R,{i,)dv 

Similarly to the fa;-case, we manipulate the equations to remove the boundary terms, 
leading to, 

+[R2.ut j dvMv.A{^^[LrR2]- + [ihm]- + [hV 
+[R2..a]- j dvMvlA{^^[LrR2]- + [th{R2)]- + [gi]-\ 

It follows that for 7^ 

j dt I ^24(0 l'< C jdt(^\ i?24(0) P +V I 4o&) I' +^ I R2vM I' II PR2 l^s 

+ II (/ - p)i?2 IlL +\ II - ^j)^2 IIL + II r IIL + II j^'hi IIL )• (4.19) 



12,. 

c / 

For ^2 = multiplying (14. lip with VzA and arguing similarly to the proof of (14.181) . gives 

j |^24(0)|'rft < C j dt(r] II PR2 III2 + II '^^l/ - P)R2 ||2,2 

+^ II 'yHi-Pj)R2 III2 + II r Il2,^ + II Wh )• (4.20) 

The only moments still to be estimated are the ?/'o-Kioiiients for .^^ 7^ 0. With C2 = 
{v'^,ip4)iv'^A,ip4)~^ and C3 = (ff , l)(f^, 1)~\ proceed similarly to the corresponding ip4- 
case discussed earlier, but start from Vz — C2VzA instead of VzA, and — c^vj instead of 
VzA. That gives 



J dt I Rl^i^z) I <C j dt(^\\ PjR, 11^2 +r] II Pi?2 III2 + II 



12.2 



+ II (/-P)i?2 11^2+^ II iyHl-Pj)R2 III2+ II r lll^ ). (4.21) 



The lemma is proved by collecting the estimates above. □ 

By using this 1-d analysis, it follows in the 2-d case that 
Lemma 4.5. The solution R2 of ^4- satisfies 



PjR2 ||2,2,2— C(^ II / ||2,2~ + II PjRl 112,2,2). 
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Proof of Lemma 14.51 We can apply Lemma|13]to -R2(0, ^z, = j dxR2{x, z, v) and, taking 
into account that gi is of order 6, get a bound for the Fourier components PjR2{0,^z), for 
6 small. As discussed at the beginning of the proof of Lemma 14.41 the components with C, 
large are under control, since the r-terms are small after division by The remaining 
components in the case 7^ 0, ^x,^z finite, are estimated by analyzing the equations for 
the moments of R2 and applying in a suitable way Lemma 14.41 The proof follows closely 
the one of Lemma 14. 4^ so we will not give all details but only point out the differences. 

Consider first the Va;-moment for ^2 = 0. Multiply the (spatial) Fourier version of (14. 7p 
by Mdv (resp. v^Mdv) and integrate. Combining the results 

£^(^2o(ex, o)^L.(ex, 0)) = ^/iex^2..(e., o)/?^.^, o) + ^/xe^^L^fe, o)^2o(ex, o) 

+i?;,^(e.,0) j dvMvzr + R2o{UO) j dvMv^VzT*. 
This gives 

j \R2vf{ix,Q)dt <C j c/t(|/?20|'(ex.,0)+r/ II Pi?2 ||2,2 + lk^(^-^)^2 ||2,2 

+l||z.i(/-Pj)i?2|||2 + lir IIL)- (4.22) 

To bound the u^-moments for 7^ 0, we use a variant of the proof in Lemma l474l Multiply 
the Fourier version of (14.71) by {{v1,vlB) — vlB)Mdv and integrate. That removes the 
hydrodynamic ^^.-term. To remove the f^r-term we also subtract the same expression 
for ^2 = multiplied with (-1)^^ We use the notation [ff := /(O - (-l)^^/(^x, 0). 
Proceeding as in Lemma [4.41 gives 

e^^{\R2o{vlvlB) - R,.B)nRl.X) = iUR2.AvlvlB) - ^2.35)^]^ 



+^^^ix([R,M. vlB) - R^^^.j^f + [i?L..J'j [^LJ' + e.i?L|(0)[^2o(t^^, vlB) - R.,.^f 
-eG[Rl,f J vlB^{M[R2r)dv + eG[R2o{i)-KiBr j Vz^{M[Rlfdv) 
+ [R;r j mQ-vIB [L7R2]" - vlB [Ih{R2)r) dv. 
It follows after integration with respect to t that for ^2 7^ 

j \R2vfmt <C j dt(\R2,f{^x,Q) + \R2,W0+V\\ PR2 ||2,2 (4.23) 

+ II v"^I-P)R2 \\l2+j2 II ^^<I-Pj)R2 112,2)- 

The moment R2vA^xi^) is under control by (I4.22p . and so the previous approach for the 
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12,2 



Waj-moment when ^2 7^ 0, gives 

J \R2vfiOdt <C J dt(\R20\\O + |^24|'(0 +^ II 

+ II uHl - P)R2 +^ II uHl - Pj)R, + II r IIL )• (4.24) 

The R2vy{C,)-'m.o'ment for 7^ is similarly treated, starting from 
d 

- (-i)^^4.,..(ex,o))(^2.,.|(0 -^2.,.|(ex,o)}. 

Inserting the corresponding right-hand sides and estimating the upcoming moments, gives 
for ^ 

j \R2vM'dt <C j dt(^\ R2.^{0) p +v II PR2 \\l2 + II iyHl-P)R2 11^2 



1 



1 



+ - II uHl-Pj)R2 \\i2). (4.25) 

For the R2vA^x, 0)-moment, use the procedure of Lemma with ifiCx^xR2v^{^x, ■) added 
as an inhomogeneous term. Thus 

e— I dz I dqTxR2o{^x,q) = ijJ'^x / dz j dqJ='xR2vA^x,q) + R2vA^x,0) 



— TT — 7T ■J —IT ■J — TT 



-27r I Mv,Txr{U-T^)dv. 



The equation for R2v^ is 



d ^ . r d 

^OlR^v, = ^fJ'CxR2v^vA^x, 0) + r^2{Q J ^^^^^2(^a', 0, v)dv. 

The resulting terms are of the same type we get before, except for 

dz j dqTxR2vAix,q) and / dz j dqJ^xR2o{^x,q) ■ 

-TV J —IT J —TV J —TV 

The former can be controlled using the observation || J-'xR2v^{^x, •) 111= \R2v^{^x,C 
together with the estimates (]4.22p . ( ]4.24p for the right-hand side. We conclude: 



2 

Zj I ! 



J \R2.f{^x,0)dt <cjdt(^ II R20 \\l2 +V II PR2 Wh 

+ II - P)R2 ||2,2 +^ II ^^(^ - Pj)R2 ||2,2 + II r IIL )• (4.26) 
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For the moment R2vy {C,x, 0) we follow the corresponding procedure of Lemma HiH The new 
terms which result from the x-derivative, R2v^vy{^x,0) and R2v^vyvA^x,0), are nonhydro- 
dynamic, and so the estimate (14.181) also holds for this moment. 

The '?/'4-moments for ^2 7^ are treated as in the 1-d case. The extra terms resulting from 
the x- derivative, contain at least one nonhydrodynamic factor. The resulting inequality is 
again KW\ . 

The moment -R24(^x,0) is obtained as in the 1-d case, but here with an additional term 
/ \R2o\'^{^x, 0)(it to the right in the final estimate (14.201) . Thus 



J \R24\\^x,0,t)dt <C J dt[\R2o\\^,„0,t) + V W PR2 111,2 + W I'Hl ~ P)R2 



1 2,2 



+ l^\\u-^{I-Pj)R,\\l^+e^\\f-\\iy (4.27) 

We now discuss the moment -R2o(^x,0) for 7^ 0. For e\a\ > cti and cti sufficiently large, 
consider the equation (14. 7p written in Fourier variables for the time and x- dependence. 
Introduce also the cutoff function jS as in Lemma HTTl Use the term ieaJ-'tJ-'xl3R2{o', C,x, z, v) 
to express the T/^o-moment. For this, project the equation along Vz — C2VzA, and along 
caV'o + v^B with C3 = — (f^, w^5) > to remove a J-'tJ-'xP R2v^-^oment. That leads to an 
expression for 

d 

TtTxl3{-iea{R2^^ - C2R2v,a) - ifJ'^xCi + ^(^20 + C2)), 

and for 

d 3 

Vx(i{ - iea{c^R2Q + R2vIb) - if^^xCs + '0^(2'^^^'^^- + '^4)), 

with (j, j = 1, ...,4 certain nonhydrodynamic moments of R2- Thus we get an expression 
for 

-leacs (^^\TtJ'xPR2o\^ + |^t^./3i?2.J') + ^ ((^t^x/3(i?2o + C2))(^C3JS:^./?(i?2.. + C4))*) . 

After division by a and integration, the boundary term is multiplied by the small coefficient 
erf ^. This can now be estimated separately at — vr and at vr using Green's formula (14.101) . 
It follows that 

II (l-X.J^t^A(-,ex, ■) Il2,2< C^(^r'(ll :FxR2{-,^x, ■) 11^2,2 (4.28) 



1 



+ - II (/-P,-)-^xi?2(-,ex, ■) 112,2,2)+ II (/-P)^.i?2(-,e. 



e 



|2 

l2,2,2 



The case of ea small requires a different argument. Consider equation (14. 7p and its 
Fourier transform in t, x and z. We denote the total Fourier transform of a function h, 
TtTxTzh by and by the Fourier transform TtTxh. With /3 defined in Lemma [4. II put 
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i?2 := PR2 and R2' := We have 

eiaR2 + ^ivJixR2 + iVzizR2 + Vzr{-1)^'' = 

eM-^Gd,^ {MR2) + e-^LjpR2 + (3Ih{R2) + e%M ■= J^- 

ei(jR2' + fivJ^^R2' + VzdzR2" = (4.29) 

eM~^Gd^^{MR2) + e~^Lj^2" + PH^2Y + eRAp' := . 

We notice that the right-hand sides contain only terms that can be estimated by contribu- 
tions either involving the nonhydrodynamic part or the hydrodynamic one multiplied by 
a small factor. 

For ^2 = we have 

eiaR2 + ^.vJ^^R2 + v^r = Af{a, 0, v). (4.30) 
We take the integral J^^ dz J'^ dq of (ICTj) 

/TT rz PTV nz 

dz I dq[eaR2'' + iiixVxR2']+iVzR2{,(^,ix,^) -'^T^ivzr{-Ti) = l dz j dqM\ 
-TT J —TV J —T: J —TV 

(4.31) 

Let M/^di'l) be a smooth function such that 

p\W'{p)fM{p)dp < 00, / p^W{p)M{p)dp = 1, / p^W{p)M{p)dp = 0, 



pW(p)M(p)rfp = 3, / pW(p)M(p)c/p = Q(.4.32) 
Jo 

Then, multiply fl4.3ip by VxVzMW{\v\) and integrate over v. The first term does not 

contribute to the hydrodynamic part of The contribution to the hydrodynamic part 
from the second and third terms are respectively 



dz j dq{tl,:,,R2')pix j dvvlvlM{\v\)W{\v\), 2m{^uR2) j dvvy,M{\v\)W{\v\). 
Let us use the polar coordinates to compute the f-integrals 



dvvtviM{\v\)W{\v\)= duoooiui dpp'M{p)W{p)=Q. 

The vanishing of the integral is due to the second condition in (14.321) . And so on the 
left-hand side there are only boundary and nonhydrodynamic terms. We put the latter on 
the right-hand side, denoted now by A/'f , 

dvVxvlM{\v\)W{\v\)^~R2{-'K) =M'. 

Vz<0 

52 



The f-integral of the ingoing part of r times fa;f^M(|f |) is zero because of the bound- 

ary conditions. In this way we have reached a control of the boundary term 7^-R2(— tt). 
Now we reproduce this term by multiplying fl4.30p by VxVzMW and integrating over Vz < 0, 

z£o-i?2..((T,e.,0)ci + z^e,i?2o(^,ex,0)ci + / dvv^vlM {\v\)W {\v\)-i- R2{-Tl) = Ml, 

Jvz<0 

where A/i incorporates J\f and all the other nonhydro dynamic terms and 

p poo P 

Cl = / dujuluz / dpp^M{p)W{p) = / dwiuliUz, 
Jsz Jo Jsz 



where S"^ = {c<j G S2,uJz < 0}, because of the first condition in (14.321) . We have used the 
second condition in (I4.32p to cancel the remaining hydrodynamic terms and as before the 
ingoing part does not contribute. The conclusion is 

Cl [eaJ2v^ [a, ^x, 0) + p^.R2o{(T, ^x, 0)] = M - ATf . (4.33) 

We will get now a second equation involving R2v^ and R20 with a similar procedure. The 
two equations together will give us the wanted estimate for these terms. 

Multiply (14.311) by f^f^Vr(|f |) and integrate over v. The first and second term do not 
contribute to the hydrodynamic part of R2. The contribution to the hydrodynamic part 
of the third term is 

t{^4 + ^Po,R2) J dvvX{tlj^ + tlJo)M{\v\)W{\v\). 

The f-integrals in polar coordinates vanish because of the second and fourth conditions in 
(I4.32p . The boundary integral is then given as before in terms of a r.h.s., denoted by N'2y 
involving non-hydro dynamic terms and A/"^, 



dvvlvlM{\v\)W{\v\)R2i-n) = Af^ (4.34) 



vz<0 



The ingoing part vanishes this time because of the second condition in (I4.32p . Now multiply 
(I4.30p by VxVzMW and integrate over Vz < 0, 

P^^xR2vM,^x,0)c2 + eaR2oi(J,^,,0)ci+ I vy,M{\v\)Wi\v\)R2i-7r)dv = X2, 

Jvz<0 

where ^ 

C2 = / duouj^uOz / dpp'^ M{p)W{p) = 3 / duu^uj 



So 



because of the first and third conditions in (14.320 . This time the ingoing part 

t I dvvlvlM{\v\)W{\v\)r,M 
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gives a contribute of order e due to the presence of the Maxwelhan M+ in the boundary 
conditions. We put this term and all the other non-hydro dynamic terms in the r.h.s. term 
denoted by A/2. By using (I4.34p we finally get 

C2eaR2nA(^, i„ 0) + c^f^^.R2o{a, ^x, 0) = ^f2 - A/?- (4.35) 

Equations fl4.33p and (14.351) give a system of two equations that allows to express R2v^ and 
R20 in terms of quantities under control provided that 02- This can be easily checked 
by direct computation. 

As an end result, it holds for ea < ai that 

II X<xii?20(-,e.,0) ||^< c(l II {I-Pj)R2 \\l2,2 + II iI-P)R2 11^2,2 +V II ^2 11^2,2 )• (4-36) 

The same estimate also holds for the i/jQ-moment when 7^ 0. The proof is simpler. 
The boundary term is this time removed by subtracting the same equation for = 
(times (— l)^"") and then multiplying first by VxVzMW and then by vIvzMW . This gives 
two equations 

e(TCi[R2vJ\- + Ci/i^^[^20]- = eaC2[R2vJ\~ + ClAi^:r[^2o]- = 5i, 

where [/]_ = /(a, e., e., ^^) - (-l)«^/(a, e., 0, i;). 

The lemma follows by collecting the previous estimates. □ 

The above study of R2 leads to 
Lemma 4.6. Any solution R2 to the problem i4.1\} satisfies the a priori estimates 

II u^I-Pj)R2 \\l2,2 <c{e II Ro III2 

+e II - Pj)g +^ II Pjg Il2,2,2 +^ II ^ Il2,2, 

II PjR2 ||2,2,2 < c(-(|| Ro ||2,2 + II ^"^(^ " Pj)9 Il2,2,2) + \ II Pj9 I|2,2,2 

+4 H 



1 



^4 II r 112, 2,~ J y 

1 1 - 

^^'R2 ||L,oo,2 < C(^ II R2 ||L,2,2 + II 7"^1 llL,2,~ +^ II llL,2,~) 

< II i?0 ||2,2 +^ II '^"^(^ - ^j)^ Il2,2,2 +JE II Il2,2,2 

+ ^ II Il2,2,~ + II ^0 ||L,2 II llL,oo,2 +^ II ^ IIL,2, 



Theorem 4.7. There exists a solution R to the rest term problem (Jjjj such that 

\R{t,x,z,v)\^M{v)dtdxdzdv < ce\ (4.37) 



[— 7r,7r] J[—7T,7r] 
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Proof of Theorem 14.71 Take the asymptotic expansion of fifth order in e. We shall prove 
that R can be obtained as the limit of an approximating sequence, and that R satisfies 
(I4.37p . Since i? is a solution to the initial boundary value problem for the rescaled rest 
term, fl4.37p in turn implies the L^-convergence to zero oi R{t), when time tends to infinity. 
Let the approximating sequence {R"'} be defined by i?° = 0, and 

+ -v^ ■ Vi?"+^ - GM-^^^^ ^ = ■^LjR^+^ + + -J(R", i?") + A, 

ot e ovz e e 

R^^^{0,x,z,v) = Ro{x,z,v), 

R^+\t,x,TT^,v) = ^ I {R''^\t,x,TTr,w) + -ij{t,x,TT^,w))\w,\Mdw 



— ^/^(t, X, =F7r, f ), xG[— 7r,7r], t > 0, ^ 0. 

Here the initial value Rq is of £-order four, A has been chosen so that j j J A{-, x, z, v)Mdxdzdv 
0, (/ — Pj)g = e{I — Pj)A is of order four, and Pjg = ePjA is of order five. 
The function R^ is solution to 

R^{0,x, z,v) = Ro{x,z,v), 

R\t,x,TT^,v) = {M-^M^){v) I {R\t,x,TTT,w) + -'ip{t,x,TTr,w))w,Mdw 

— ipitjXj^TijV), xE[—n,7i], t > 0, 

Split R^ into two parts Ri and R2, solutions of (14.41) and (14.71) . respectively, with g = eA. 
Then using the corresponding a priori estimates. Lemma IT3] and Lemma HSl together with 
the exponential decrease of ip, and the e-orders of -Rq and A, we get for some constant Ci. 

II I^^R^ ||oo,oo,2< CiE^, II ^^R^ l|2,2,2< CiE^ , 

By induction 

II J^^R^ ||oo,oo,2 < Sci^t, j <n + l, 

II Z/^(i?"+l - i?") ||2,2,2 < 026^ II J^HR'' - ^""') l|2,2,2, n > 0, 

for some constant C2. Namely, if this holds up to n*^ order, then 

ot e ovz 

e-^ e e 

{R''+^ -R''+^){0,x,z,v) = 0, 

(i?"+2-i?"+i)(t,a;,T7r,t;) =M-iM^ [ - i?"+i)(t, a;, ^vr, M;)|w;,|Mrf«; 

X G [— TT, tt], t > 0, Vz ^ 0. 
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Here 



It follows that 

< C£-5 ( II Z/5i?"+l ||^^^^2 + II ^^^i?" ||oo,oo,2 ) || - R^) 

Consequently, 

II ||2,2,2 <|| - ||2,2,2 +...+ II //^(i?' 

for £ small enough. Similarly || K^^'^ || 00,00, 2< 2ci£5. In particular {K^} is a Cauchy 
sequence in L|f([0, +0x3) x Q x R^). The existence of a solution i? to (14. ip follows, and 
the estimate (14.371) holds. This means that there is a sequence of Lebesgue points in time, 
{tj}^i with tj tending to infinity with j, where the || ■ ||2,2- norm of the solution R tends 
to zero. But the ||-R(t, ■ )||2,2 for fixed t > tj is uniformly bounded by the norm at tj plus 
some tail integrals from tj to 00, hence tends to zero when time tends to infinity. □ 
This completes the study of the i?-term and the stability theorem follows. 
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